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Duct  chord 

Mean  line  ordinate  of  the  duct  section  measured  from 
the  nose>taiI  line 

Complete  elliptic  integral  of  the  second  kind 
(o/2  chord-diameter  ratio  of  the  duct 
Complete  elliptic  integral  of  the  Hrst  kind 
Modulus  of  the  elliptic  integrals 
Ring-source  strength 
Duct  radius 

Cylindrical  coordinates 

Half-thickness  ordinate  of  the  duct  section 

Free-stream  velocity 

Axial  component  of  induced  velocity 

Component  of  freo-stream  velocity  in  direction  of  duct 
axis 

Radial  component  of  induced  velocity 

Radial  coordinate  nondimensionalired  by  the  propeller 
radius 

Axial  coordinate  nondimensionalized  by  the  duct  chord 

Angle  of  attack  of  a  duct  section 

Ideal  angle  of  attack  of  a  duct  section 

Relative  angle  between  freo-stream  velocity  and  duct 

Ring-vortex  strength 

Mass  density  of  fluid 

Angular  velocity 
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Subscripts 


d 
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fL 

d<}> 


Duct 

King  source 
Ring  vortex 

Trailing  vortex  system  of  the  vortex  cylinder 


Note:  Many  functions  are  defined  in  the  text. 
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ABSTRACT 


A  coinputor  program  is  prt'sontpd  which  calculates  the  aerodynamic  char¬ 
acteristics  of  annular  airfoils  on  an  lBM-7()90  high-speed  computer.  A  brief  re¬ 
view  of  the  theory  is  also  presented.  Experimental  and  computer  results  indi¬ 
cate  that  tile  theory  gives  reasonable  prediction  of  the  lift,  induced-drug,  and 
moment  coefficients  and  also  of  the  pressure  distribution  except  when  separa¬ 
tion  is  presort  on  the  annular  airfoil.  The  computer  program  can  also  be  used 
for  the  desigi  of  ducted  propellers  if  an  infinite  number  of  blades  is  assumed. 

ADMiNiSTRATtVE  INrORMATiON 

This  work  was  covered  by  Subproject  S*R011  01  01  of  Task  0401  under 
the  Bureau  ot  Ships  In-Ilousc  Independent  Research  Program. 


1.  INTRODUCTION 

Annular  airfoils  have  found  use  in  both  aeronautics  and  naval  architecture  where 
their  application  has  generally  been  as  shroud  rings  around  propellers.  The  aeronautical 
application  has  been  t)  increase  the  thrust  of  propellers  for  hover ‘..ig  flight  such  as  applied 
to  vortical  *ake-off  craft.  In  naval  architec-urs  Kort  nozzles  have  been  used  for  many  years 
for  increasing  the  efficiency  of  heavily  loaded  propellers,  whereas  pumpjets  are  a  more  re¬ 
cent  innovation  for  delaying  cavitation. 

in  mo.st  applications  the  annular  airfoil  is  used  as  an  integral  part  of  the  propuisor, 
and  a  theoretical  treatment  must  con.nder  the  interaction  between  the  propeller  and  airfoil. 
The  present  theories,  however,  use  an  interative  procedure  between  propeller  and  annular  air¬ 
foil  theory  so  that  the  usefulne.ss  of  the  complete  theory  d.-pends  on  the  verification  of  each 
part.  One  purpose  of  this  report  is  to  present  a  comparison  of  theoretical  and  experimental 
results  (or  a  number  of  annular  airfoils. 

The  nr.sl  theoretical  development  of  annular  airfoils  is  that  due  to  Dickmann'  wlo 
represented  the  annular  airfoil  by  a  distribution  of  ring  vortices.  This  approach  implies  that 
the  airfoil  has  no  thickness.  Other  contributors  have  been  Stewart.^  Kuchemann,^*'*  Pivko,'*’ 
Malavard,^  and  Hacques.^  A  more  complete  review  of  each  of  these  contributions  can  be 
found  in  References  8  and  9.  The  most  complete  theoretic.al  approach  to  date  has,  however, 
been  that  by  Weissinger*®' *  *  who  extended  the  work  of  Dickmann  and  has  considered  ducts 
with  thickness,  of  arbitrary  shape  and  at  an  angle  of  attack.  His  mathematical  model  of  the 
annular  airfoil  was  a  distribution  of  ring  vortices  and  ring  sources  lying  on  a  cylinder  of 
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Reference*  are  listed  on  page  82. 


diameter  representative  of  the  duct  and  of  length  equal  to  the  ducL  This  means  that  the 
boundary  conditions  are  linearized  and  are  satisfied  on  a  representative  duct  cylinder  and 
pot  on  the  duct  surface. 

A  number  of  investigators  have  used  the  Dickmann-Weissingcr  tnathematical  model  a 
have  developed  numerical  methods  for  computing  the  pressure  distribution  on  annular  airfoi 
Bagley,  Kirby,  and  Marcer*^  have  presented  a  simplified  method  where  the  boundary  condi* 
tions  are  satisfied  at  only  a  few  points  along  the  chord.  A  more  complete  method  is  that 
given  in  References  13  and  14,  which  has  presented  quite  an  extensive  set  of  tables  for  a 
limited  r.u  nber  of  airfoil  configurations. 

In  addition  to  the  comparison  of  theoretical  and  experimental  results  mentioned  pre¬ 
viously,  this  report  presents  a  computer  program  for  the  IBM*7C90  which  enables  not  only  t 
pressure  disuibution  of  the  annular  airfoil  to  be  obtained  but  also  other  aerodynamic  chara* 
teristics.  Specifically,  the  program  calculates  the  pressure  distribution,  both  line.ar  and  nc 
linear*^  for  an  annular  airfoil  of  arbitrary  section,  the  ideal  angle  of  attack,  the  velocity  di 
tribudon,  both  on  and  in  the  flow  field,  and  the  various  aerodynu.’nic  forces  which  arise.  T 
annular  airfoil  is  assumed  to  be  axisymntetric  but  nay  have  an  angU-  attack  and  be  in  an 
axisymmetric  flow  field.  Considerati.in  of  such  a  flow  field  allows  the  inclusion  of  an  infi- 
nitely  bladcd  propeller  or  a  central  b>dy  within  the  ducL  The  mathematical  model  used  for 
these  calculations  is  the  Dickmann-Heissinger  model  and  the  approach  is  that  of  Reference 
9  and  15. 

The  following  discussion  is  divided  into  four  main  sections.  The  development  of  th 
theory  is  reviewed  briefly,  then  the  computer  program  is  described,  and  finally,  the  results 
of  the  calculations  and  comparison  with  the  experimental  results  are  presented. 

2.  LINEARIZED  THEORY  OF  THE  ANNULAR  AIRFOIL 

The  linearized  theory  of  the  annuiar  airfcil  has  been  developed  adequately  in  the  rel 
erences  cited,  but  a  brief  development  will  be  repeated  here  for  completeness.  Details  of  l 
mathematical  manipulations  which  have  been  omitted  will  be  found  in  References  9  and  15. 

As  indicated  in  the  Introduction,  the  method  of  singularities  is  used  for  ropresenuti 
of  the  flow  field  about  the  annular  airfoil.  The  mathematical  model  used  will  be  a  distribu¬ 
tion  of  ring  vortices  and  ring  sources  lying  on  a  cylinder  of  a  diameter  representative  of  ih« 
duct  diameter  and  of  length  equal  to  the  duct  length  (the  Dickmann-Weissingcr  model).  In  t 
use  of  this  approach  a  number  of  assumptions  are  implied.  Briefly  those  are: 

fu  The  fluid  is  inviscid  and  incompressible  and  no  separation  occurs  on  the  duct. 

b.  B>>dy  forces  .such  as  gravity  are  neglected. 

c.  The  free-stream  flow  is,  in  general,  axisymmetric  but  may  have  a  small  cross- How 
component. 

d.  The  annular  airfoil  is  axisymmetric  and  of  finite  length. 
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o.  The  annular  airfoil  can  be  represented  mathematically  by  a  distribution  of  ring  vortices 
and  ring  sources  along  a  cylinder  of  constant  diameter.  This  implies  that  the  boundary  condi¬ 
tions  are  linearized. 

f.  The  trailing  vortex  system  of  the  annular  air.'cJl  has  the  constant  diameter  of  (he  an¬ 
nular  airfoil  and  extends  from  the  annular  airfoil  to  infinity. 


2.1  BOUNDARY  CONDITIONS 

The  coordinate  system  used  is  a  cylindric.!!  system  (r,  Z)  with  the  axis  located  at 
the  trailing  edge  and  on  the  centerline  cf  the  duct.  For  ronvenience,.the  axial  coordinate 
will  be  nondimensionalized  by  the  chord  c  and  the  radial  coordinate  by  the  reference  radius 
of  the  duct  Rj.  Figure  1  shows  the  annular  airfoil  coordinate  system  and  Figure  2  is  a  de¬ 
lineation  of  the  system. 

The  cro.«s  section  of  the  annular  airfoil  is  assumed  known  and  will  be  delineated  by 
a  thickness  distribution,  a  camber  distribution,  and  an  angle  of  attack. 

In  terms  of  the  thickness,  camber,  and  angle  of  attack,  the  slope  of  the  outer  surface 
of  the  duct  is  given  by 


u'(Z)  =  C'(Z)  *  tan  «  +  S'{Z) 

and  the  slope  of  the  inner  surface  by 

6  (Z)  =  Cj'{Z)  ♦  tan  a  -  S\Z) 


The  boundary  conditions  to  bo  satisfiisi  oi.  the  surface  of  an  annular  airfoil  are  that 
the  normal  velocity  must  be  zero  and  that  the  Kutta  condition  must  be  satisfied  at  the  trailing 
edge.  In  linearized  theory.  Reference  9,  this  means  that  the  radial  velocity  at  'he  reference 
cylinder  must  be  c«{ual  to  the  slope  of  the  section,  or 


r; 

y  {Xj  i  0,  Z)  - -iq(Z)  *  tana  i  S'(2)l  (0  ^  Z  ^  1)  [2.i.lj 


And  at  the  trailing  edge  the  radial  velocity  is  zero,  or 


(A;,i0,  o,  0)  =  0 


(2.1.21 
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2.2  DERIVATION  OF  THE  RING  VORTEX  AND  SOURCE  DISTRIBUTIONS  TO  REPRESENT 
THE  ANNULAR  AIRFOIL 


In  order  to  calculate  the  pressure  distribution  on  the  annular  airfoil,  the  flow  field  a- 
round  the  foil,  and  the  aerodynamic  forces,  it  is  necessary  to  obtain  the  strenath  of  the  ring 
vortices  and  ring  sources.  The  ring*vortex  strength  and  ring'Source  strength  must  he  of  such 
magnitude  that  they  induce  radial  velocities  which  satisfy  the  boundary  conditions  given  by 
Equation  [2.1.1].  By  substituting  the  equations  for  the  radial  velocity  induced  by  the  ring  vor> 
tices  and  sources  into  Equation  [2.1.1],  a  singular  integral  equation  is  obtained  which  can  be 
solved  for  the  various  singularity  distributions  as  shown  in  the  following  discussion. 

The  nondimensional  elementary  circulation  of  the  ring  vort.ices  v,;ll  be  taken  to  be 
y(^,  Z)  and  the  strength  of  the  ring  sources  to  be  ^(<5.  Z).  If  the  ring-vortex  strength  is  a 
function  cf  t'-  <  angular  coordinate,  then  from  vortex  theory  a  trailing- vortex  system  exists 

1  dy 

behind  the  duct.  The  strength  of  the  elementarv’  vortex  is - ,  and  this  vortex  system 

induces  radial  velocities  at  the  duct.  Since  the  velocities  are  linear,  they  are  additive,  and 
the  radial  velocities  by  the  duct  are  found  to  be 


(O^Z^l)  [2.2.1] 

.And  the  various  terms  are  defined  as  follows: 

The  radial  velocity  induced  on  the  duct 
by  the  ring- vortex  system. 

The  radial  velocity  induced  on  the  duct 
by  the  rinj.^  source  system. 

The  radial  vv-’oeity  induced  on  the  duct 
by  the  trailing-vortex  system. 


[»;(A^,  si,  Z)iy 

Z)]y 

z)]^y 


cos 


An  arbitrary  axisymmetric  vplof  ity  at  the 
duct  such  as  induced  by  an  infinitely 
bladed  propeller  or  a  center  body. 

The  radial  component  of  the  frec*stream 
velocity  when  the  duct  is  at  an  angle  of 
attack. 


The  radial  velocity  components  induced  by  the  various  singularity  distributions  at  the 
reference  cylinder  have  been  derived  in  Reference  9  and  are 


For  the  ring  vortices: 


pi 

1  ,12.^, 

y-ix,,  Z] 

-*  y  ~  0  0  [4/(^(Z-Z'l^  +  4  sin^ 

[2.2.21 

For  the  ring  sources: 


1  7n 


kff 


(l-cos(<i-ci')l  g(6:z‘) 


0  0  *  4  .sin2_(p-^'))3/2 


Z) 

19.2.31 


For  the  trailing- vortex  cylinder: 


cot-;j-(si-«6') 


\^Hz-z' 


2f4,Z~Z') 


1 


')2.4sin2~(0.6') 

•7 


12.2.41 


w  here  ft 


'2R. 


These  radial  \clocities  are  substituted  into  ihe  boun  !ary  condition  12.1.11  which  leads 
immediately  to 


^6.  Z)  -  -2  5'(Z)  12.2.51 

From  this  equation  it  is  seen  that  the  strength  of  the  source  rings  is  a  function  only  of  the 
thicknes.s  slope  and.  for  an  axi.symmetric  duct,  is  independent  of  angle.  In  linearized  airfoil 
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theory  the  same  equation  is  obtained.  Substituting  the  slope  for  the  source  strength  into 
Equation  (2.2.3],  the  radial  velocity  induced  by  the  ring  sources  at  the  cylinder  becomes 


^  0,  Z)  -  -iy*  I  AW  -  £:(^)]  dZ ' "  5  '(Z)  12.2.6] 


.’here  =- 


1 


and  A^A)  and  £(A)  arc  complete  elliptic  integrals  of  the  fir.st  and 


h\Z-Zy  M 

second  kind,  respectively. 

Substitution  of  the  various  radial  velocities  into  the  boundary  condition.  Equation 
[2.2.1],  led  to  the  ring-source  strength  given  by  Equation  [2.2..'i].  This  same  substitution 
also  led  to  the  following  integral  equation  for  the  ring*vurte,\  circulation: 


1  2;r 


2A<Z-Z')  cos(^-<5')  >'(P,'Z') 


‘d6‘dZ' 


0  0  [4A2(Z-Z')2  ♦  4  sin2— (,;S-cS')) 

2 


3/2 


.  1  2ff  ^ 

ACC  1  2A(Z-Z') 

dJ  ' rr=r=Tf — 

0  0  ^'4A2(Z-Z')^-r4  sin^— (sS- cS') 


dy 

—  (6:Z')dif>'dZ' 
d6 


f  hr  K  1 

-  jlC;(Z)  *  tan  o  4  «,  cos  <i] - j  S'(Z ')  A[H{A)  -  E{A)]  dZ  -  d  z)\ 


-  f/(sS,  Z)  12.2.7) 

This  equation  is  a  singular'integral  equation  for  the  ring>vortex  circulation  y(d-.  2). 
In  order  to  solve  this  equation,  an  as.sumption  is  made  th.at  y(0,'  Z ')  can  be  expander,  in  a 
Fourier  series  in  6,  i.e., 


y(6,'  -  ')  «  y,(Z ')  cos  n6  *2’  A^(Z ')  sin  n«^'  (2.2.81 


n  =  0 


n  >  1 


6 


The  function  6'(^,  Z)  is  also  expanded  in  a  Fourier  series  in  ifi,  •*-  is  immediately  obvi 
ous  because  of  the  form  of  f/(^,  Z}  thut  only  two  terms  of  the  series  exist,  i.e., 

Z)  -  -i  cos  12.2.9) 


where 


f  A  /•*  1 

«o  =•  A\C;iZ)  >  tan«]  — -y  S'(Z')  k[K{lc)-E(k)\  a?  {X^,  Z)| 

* 


If  Equations  {2.2.F0  and  [2.2.9]  are  substituted  into  Equation  [2.2.7]  and  the  neces¬ 
sary  icnthemat.cal  manipulations  are  performed,  a  singular  intcgrrl  equation  is  obtained  for 
the  Euler  coefficients  ff^iZ)  arid  ff^(Z).  All  other  Euler  coefficients  are  zero  since  a  Fourier 
series  is  unique  and,  hence,  the  coefficients  can  be  equated.  Only  and  u^exist  on  the 
right-hand  side  of  Equation  [2.2.7].  The  two  integral  equations  obtained  are 


/ 


a(Z-Z') - f!(Z) 

UZ-Z) 


[2.2.10] 


f  - 7-  Hj(Z -2  ')dZ  '^nh  /  g^{Z  ')dZ'^  -2n«,  [2.2.11] 

d  (Z  — Z  )  d 


where 


[2.2.12] 


y(Z-Z')  =  Jt|4A2(Z-Zy  [A(*)-E(*)]  -  2E(*>| 

W  ‘ 

H{Z)  =4a[<7,'(Z)  ian«l  -4.r  Z)  -4A  J  S'{Z')HK(k)-E(k)]dZ'  [2.2.13] 


A 


1 

h\z~zy  ♦  I 


7 


rCZ-Z')  -i-U2-X*2)2  E(k)  -4(1-X2)2 

Solution  of  these  two  integral  equations  gives  the  circulation  distribution  for  an  axisymmctric 
annular  airfoil  at  an  angle  of  attack,  i.e., 

Z)  =  g^JiZ)  g^Z)  cos  ^  (2.2.141 

Actually,  soluusn  of  Equation  (2.2.10}  gives  the  circulation  di.stribution  for  an  axisymmctric 
annular  airfoil  at  zeic  an^te  of  attack.  It  should  be  noted  that  the  thickne.°^'<  distribution  and 
the  arbitrary  axisymmctric  radial  velocity  each  contribute  to  the  strength  of  the  ring  vortex. 
Solution  of  Equation  (2.2.11]  gives  the  circulation  distribution  on  a  right  circular  cylinder  at 
an  angle  of  attack  with  length  equal  to  the  duct  chord  and  diameter  equal  to  the  representa¬ 
tive  diameter  of  the  duoL  This  circulation  distribution  is  independent  of  the  .shape  of  the 
duct. 

The  equations  for  the  circulation  distributio:  s  have  been  derived  in  terms  of  elliptic 
integrals.  They  can  also  be  derived  in  terms  of  the  half-oder  Legendre  function  of  the  sec¬ 
ond  kind.^  If  the  annular  airfoil  is  not  cxisymmetric  and  the  propeller  has  a  finite  number  of 
blades,  the  uso  of  the  Legendre  functions  is  preferred. 

In  the  next  few  sections  the  two  singular  integral  equations  just  derived  are  reduced 
to  Fredholm  equations  of  the  second  kind  and  the  solution  of  these  equations  is  discussed. 

2.3  REDUaiON  OF  THE  INTEGRAL  EQUATION  FOR  THE  AXISYMMETRIC  CASE 

Following  Muskhotishvili,**  as  shown  in  Reference  9,  Equation  (2.2.10)  can  be  re¬ 
duced  to  a  Fredholm  equation  of  the  second  kind  by  adding  and  subtracting  the  tenr.  g(Z-Z') 
from  the  kernel.  Since  at  Z  ■  Z'the  kernel  is  equal  to  -2,  i.e.,  g{Z-Z')  -  g{0)  =  -2,  Equa¬ 
tion  (2.2.10]  becomes 


J  I  c  1 


The  integrand  of  the  integral  on  the  right  is  not  singular  at  the  point  (Z  *  Z')  and  has  the 
value  of  «ero.  This  equation  is  in  the  form  of  the  well-known  Cauchy-type,  singular-integral 
equation  and  has  a  unique  inverse  given  by 


(2 


'±^f^(Z')dZ'*2j  g^iZldZ' 


(2.3.2] 


8 


I 

f 


To  this  equation  the  secoa'*  boundary  condition,  Equation  {2.1.21,  is  applied.  This  boundary 
condition  implies  that  the  circulation  at  the  trailine  edge  must  be  zero,  and  the  constant  term 
1 

j  g^Z.  ')dZ  '  is  chosen  so  that  ^(,(0)  =0.  If  f^iZ ')  is  substituted  inU>  this  equation  for  /^(Z ') 
0 


and  the  constant  chosen  as  just  descril^jd,  a  Fredholm  equation  of  the  second  kind  is  ob¬ 
tained  for  the  circulatior  distribution: 


1 

=•  v'  *  *  J  2 


dZ' 


[2.3.31 


where 


m 


Z' 


//(Z ')  dZ ' 


2n^  JF  Z"  (Z"-Z)fZ"-Z') 


A  new  dependent  variubie  is  defined  by  this  equation  from  which  the  circulation  di.-'tri button 
can  easily  be  obtained.  This  redefining  is  necessary  since  the  circulation  distribution  gJiZ) 
is  generally  infinite  at  the  leading  edge.  The  case  where  gQ(Z)  Is  not  infinite  at  the  lead! 
edge  will  be  discussed  later. 

Both  functions  f(Z)  and  A'j(Z,  Z')  are  Cauchy  Principal-value  integrals  and  to  evalu¬ 
ate,  part  of  the  integrand  will  be  expanded  in  ,i  Fourier  series.  For  convenience,  a  change 

of  variable  of  the  form  Z  »  — (1  *  cos  0)  is  made  in  Equation  [2.3.31.  Then 


r 

Fo*(<0  =  Fo(«)  sin  -i  t?  =  f{0)  *  J  h\0.  O')  g;(e')  dd'  [2.3.41 


9 


where 


m 


1  1  / 

- cos  —  0  C 

2„2  2  / 


(l-co»  <?') 
(cos  O'- cos  0) 


H{e')dO’ 


(2.3.5) 


11  1  / 
K{0,  0 ')  =  —  cos  —  0  cos  — 0 '/ 
-2  2  2  J 


,/  (1-cosO") 

2  +  gico'  9"-cos0')| 

/  (cos  0 cos  0) 

0 

cosO  '-cosO'  J 

\d9"  (2.3.6) 


To  evaluate  the  function,  f(0),  //(O')  is  obtained  in  series  form.  This  function  is  ^iven  bv 
Equation  (2.2.13)  with  the  proper  change  in  variable.  The  following  half-range  Fourier  .serie.s 
expansions  are  made  for  the  slope  of  the  thickness  distribution  and  the  section  mean  line: 


~  2)  =  C?j'(0)  -  —  (A^,  0)  -  cos  mO  (2.3.7) 


01  >•  1 


where 


and 


S(Z)  =  S(0(Z)1  .sin  ,7i0 

mm  I 


(2.3.8) 


where 


sin  mO  dO 
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The  slope  of  ihe  thickness  distribution  is  then  given  by 


dS  2  ■  ^ 

S'(Z)  - — {e(Z)\  =  -  cos  md 

dZ  sin  " 

1 

Use  of  the  Fourier  scries  expansion  causes  no  restrictions  on  the  th. cl' ness  distribution; 
however,  Equation  [2.3.7]  requires  that  the  camber  line  slope  bo  finite  everywhere,  intro* 
ducing  the  expansions  into  Equation  [2.2.13]  yields  the  following: 


«»  m  ««  _ 

H{e‘)  -  |4a  ^un«  +  gp  cos  me  yAl^ 

*1^1  m  K  1 


[2.3.9] 


where 


G{6:m) 


n 

'*3/ 


-f  t^k  {K{k)  -  E(k)]  cos  m  (d 1^)  dt 


[2.3.10] 


*2 


4 

A2[cOS  0'-COS(d'-f*})2  ♦  4 


The  elliptic  integral  of  the  first  kind  K{k)  has  a  logarithmic  singularity  at  A  -  1,  which  re¬ 
sults  in  the  integrand  of  the  left-hand  integral  of  Equation  [2.3.9]  having  a  logarithmic  sin¬ 
gularity  at  this  point.  This  is  a  singularity  which  can  be  removed  by  making  the  change  in 
variable  of  0’  -  d"=  t\  This  has  been  done  in  the  right-hand  integral. 

To  complete  the  solution,  the  function  G{0' m)  is  expanded  in  a  Fourier  cosine  series 
in  f.'i.e.. 


G(0;  ffi)  =  cAm) 


'3' 


(m)  cos  pO’ 


P-1 


(2  3.I1J 


11 


1 


where 


It 

aAm)^-f  G(0:m)d0' 

2  / 

a(m)~—  I  G(d:  m)  cos  pO'dO’ 
P  n  J 


With  this  expression  for  G(0^ m).  Equation  [2.3.9]  is  substituted  into  Equation  [2.3.3) 
and  the  integral  for  f(d)  cru)  now  no  evaluated.  After  the  order  of  integration  and  summation 
are  interchanged,  the  resulting  integrals  are  of  the  Glaucrt  type  and  f(0)  then  becomes 


f(0)  =  j^2(Un«  e  C^)  e cosjt? 

[-2^„  sin  m0  sin  jO 

L  1  r.  *.  1  J 


[2.3.121 


where 


6A  ^ 

- - m  ^a  (m)  sin  p0 

P-i 


To  obtain  the  kernel  K{0,  O')  also  involves  the  evaluation  of  a  Cauchy  principal-value 
integral.  The  method  used  to  evaluate  this  inu-gral  is  to  expand  the  part  of  the  integrand 


[2  +  ^(cos  (?"-ro3  C')  1 
cos  9 "-cos  0'  J 


in  a  half-range  Fourier  cosine  series  in  01  This  term  is  continuous  everywhere  for  0-0" 
and  has  the  value  zero  for  0.' 


i 

t 


where 


i  +  e(cos  (?''-cos  0')1  '177 

— — - r-^  -  +  yf>AO')  cos  nd" 

cosfl“-cos«'  J  ° 

«  •  I 


12.3.131 


1  f  fo -f  c(cos  i>"-cos  W')! 

^ - :  do" 

®  r  J  L  cos  0"-cos  0'  J 
0 

.  2  t  r2  +  ff(cos  0"-~cos  0')"1 

b(e')^—l  — — - -\  cos  nO" do'- 

"  L  0'‘-cos  0’  J 


Substituting  this  equation  into  Equation  f2.3.6)  and  evaluating  the  integrals  gives  the 
following  for  the  kernel  K{0,  0')\ 


K(0.  O’)  -  —  cos  ^0'  -i^iO’)  cos-i  d  +  sin  -^0^(0')  pin  n0  12.3.141 

*-  (1^  1 

Both  f(0)  and  KiO^O')  are  now  in  a  form  which  can  bo  handled  numerically  so  that  the  Fredholm 
equation  of  the  second  kind.  Equation  (2.3.3],  can  be  solved  for  the  circulation  distribution. 
Several  method.s  exist  for  the  solution  of  this  typo  of  integral  equation;*^  however,  the  ker* 
nel  A'(d,  O’)  is  of  the  degenerate  (or  product)  type  and  the  method  applicable  to  this  type  of 
kernel  will  be  used  here.  Details  of  this  method  are  given  in  Reference  17  and  only  resulting 
equations  arc  given  here.  Following  this  procedure,  the  equation  for  the  circulation  di^tribu* 
tion  becomes^ 


=  /(0)  +  C(0)  /?g(0)/lQ  +  C(d)^jSin  0  +  . . .  +  C{{?)/l^sin  nO 


12.3.151 


where  /(d)  is  given  by  Equation  (2..3.12] 


and  is  given  by  the  following  set  of  simultaneous  equations: 

“^^01  “^2^02" - “4) 

*•^0^10  +4[('“^u)  "■^2^I2~ - 


12.3.16] 


whore 


‘I 


U 

'  Tnf  /  -  0,  1,  2,  . . .  n) 


DJd’)  m  -cot  -  d' 


/><«')  .» sin  fO' 


acd 


•  —  -1  *0  — 

5.  !>(Un..C„)^F,  (■•-O.I.2...n) 

»  «•*!*  «*l  mm  I 


V 

fi^J  6,<0')  OO8^J0'd0' 


n 

d.^  m  j  b^{0')  sin  e'sin  in0' d0’ 


TT 

fim  - /  »«" 
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•V 


T 


2.4  THE  IDEAL  ANGLE  OF  ATTACK 

In  didcussing  Equation  (2.3.2]  of  the  previous  section,  the  statement  was  made  that 
the  circulation  distribution  of  the  annular  airfoil  was  generally  infinite  at  the  leading  edge. 
This  is  not  the  case,  however,  if  the  annular*airfoil  section  is  at  its  ideal  angle  of  attack 
i.e.,  if  the  stagnation  point  occurs  at  the  leading  edge.  This  angle  can  be  obtained 
from  Equation  [2.3.2]  by  not  only  iraking  y^{0)  >  0  but  also  by  taking  the  circulaticn  at  the 
leading  edge  to  be  zero,  ^^(1)  «  J,  and  solving  ^or  the  angle*: 


dZ' 


[2.4.1] 


where 


If  this  value  of  the  ideal  angle  is  substituted  in  place  of  the  section  angle*,  an  equation  is 
obtained  for  the  ideal  circulation  distribution: 

1 

Yij  -  fjZ)  *  f  K^/Z.  Z  )  yJZ')  dZ‘ 

0 


where 


[tf(Z')-4ff  tan*] 

(Z'-z)vzvr^ 


1  r -  /  2  +  p(Z"-Z') 

K  AZ,  Z')  -  —  a/Z(1-Z)  # - — - == 

2n^  J  (2"-<Z')(Z"-2) 


y/Z’ll-Z") 


rdZ" 


In  this  equation  it  is  not  necessary  to  solve  for  a  pseudo>circulation  since  y,,^(0)  not  sin¬ 
gular  at  the  lending  edge.  By  the  following  the  procedure  of  the  previous  section,  the  coef¬ 
ficients  ffjiZ)  and  K^j{Z,  Z')  are  obtained  as  follows: 


15 


I 


\ 


e')  *  —  sin  sinn0 


A  a  1 


The  coefficie?it8  given  in  thA  previous  section.  The 

ideal  circulation  distribution  is  given  by  Equation  [2.3.15]  wiih 


m-- 


Pj(fl)  .  0 


•kd  is  given  by 


“ . 

-^Cjj  +^(1-C22)  - . 


. 

All  the  coefficients  are  the  same  es  given  in  the  previous  section  e^:ept  for  d.,  which  is 
given  by  ^ 

—  m 

M  *  t.  mat 


le 


1 


f 


Once  the  ideal  circulation  distribution  is  obtained,  the  ideal  angle  of  attack  at  which  the 
section  must  operate  for  the  stagnation  point  to  be  at  the  leading  edge  is  given  by  Equa* 
ticn  (2.4.11.  Substituting  for  '••..to  this  equation  yields  *he  following  for 

•a  -  j;/  y^')  sin  (2.4.31 

2.5  REDUCTION  OF  THE  INTEGRAL  EQUATION  FOR  THE  ANNULAR  AIRFOIL  AT  AN 
ANGLE  OF  ATTACK 

In  Section  2.2  an  integral  equation  w*s  deilved  for  the  circulation  distribution  of  an 
annular  airfoil  at  an  angle  of  atta.  k.  This  equation  showed  for  the  linearized  theory  that  the 
section  shape  had  no  effect  on  the  circulation  distribution,  which  was  dependent  only  on  the 
chord-diameter  ratio  h  and  the  angle  of  attack 

Thd  rcdu^.tion  of  this  Equation  [2.2.11]  to  a  Fredholm  equation  of  the  second  kind  fol¬ 
lows  the  procedure  outlined  in  Section  2.3  and  is  more  fully  described  in  Reference  9.  The 
resulting  equation  is 


g*{Q)~  (sin^o)  y,(») 


where 


-2«,cos~0-.lJ  ^cos^d' ros 
(sin^d)^^(d^)  sin  mdjj  dd' 

ir  1  r  r d"-cos  d')  -1  I 

A+-/  - - \de- 

2  »./  I  (cos  d"-cos  O')  I 

2  /  ^(cos  d"-cos  d')  -l1 

(d')  —  /  - - - - - —  cos  md"d6" 

(cos  e  -cos  d  )  I 


(2.5.11 


I 
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Here,  as  in  Section  2.3,  it  is  necessary  to  solve  for  a  pseudo-circulation.  From  Reference  17 
the  solution  of  this  equation  for  the  pseudo-circultfton  is  obtained  in  the  following  form: 

12  12 

»  2«,  cos— - ilftCos  — sin  d+  -sin  2d+. . .  +  4. sin  ntf]  (2.5.2] 

’  '  2  n  ^  2  n  ^  2  " 

The  coefficients  are  obtained  from  the  set  of  simoltaneous  equations  given  by  Equa¬ 
tion  [2.3.1],  except  that  the  coefneients  and  m- 


V 

C..  •  0.(0')  sin  d'dd';  S  /=0,  1,  2  ...  n) 

0 

n 

d.  =  2m^J  cos^  0' bJiO')  d0'\  1,  2,  3  . . .  n)' 


-cot-fl' 
«  sin  jO'; 


(y»  1,  2,  3  . . .  n) 


3.  VELOCITY,  PRESSURE.  AND  FORCES 
3.1  VELOCITY  DISTRIBUTION 

!n  linearized  theorv'  the  flow  field  is  given  by  summing  the  free-stream  velocity  and  the 
velocities  induced  by  singularities  in  the  flow  field.  For  the  annular  airfoil  at  zero  angle  of 
attack,  this  means  that  to  the  free-stream  velocity  must  be  added  the  velocity  induced  by  the 
ring  sources  and  ring  vortices.  For  the  zero-incidence  case,  only  axial  and  rad'al  velocities 
are  induced  by  the  duct  and  these  are^ 

Axial-induced  velocity: 


k\cc3  0-COS  0’)^ 


K(k^)-E{k^) 


cos-(0^+e)e^  dP  [3.1.1] 
2  2 
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and 


2  ■  •^(1  +  cos  0) 


[3.1.5] 


The  foregoing  equations  for  the  induced  velocities  have  been  derived  in  a  form  for  ease 
in  making  numerical  calculations.  Changes  in  variables  have  been  made  in  order  u>  remove 
singularities  from  the  integrand. 

On  the  anrulf>r  airfo:l  itseif,  these  equations  are  considerably  simplified  since  X  =  A^. 
However,  across  a  singular  vortf.  sheet  there  is  a  discontinuity  in  the  tangential  velocity 
while  across  a  singular  source  sheet  there  is  a  discontinuity  in  the  normal  velocity.  On  the 
duct  itself  then;  i.e.,  A'  «  Xj  and  0  -  Z  -  1, 

Axial-induced  velocity: 


X  1  h  1 

-  -  j  k[K(k)^E(k)]  cos  ^  (0^-0)  A6^d3  “  [3.1.6] 

-  -Jv  ^  •'iTTr-  ^ 


and 


•y  -inr 


^)]  -  -  f  ^  p 

I  ^  J  nj  (cos  0-cos  0')  sin 

V  •«  n  =  1 


.sinn0  13. 1.7] 


n 


(0')  kE{k)  cos  n6'  do ' 


Radial-induced  velocity: 


Pr  1  1  /  SS^n'k  /, 

I  xMd’  I  “  '~T  / - a - 0-COS  0 

L»  njJ  (cos  0-co.s  0  ) 


')2  lE(k)-E(k)J-E(k) 


^cos  — 0'rf0' 
o 


—X 

sin  Oja^ 

« •»  1 


sin  n0 


[3.1.8] 
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and 


13.1.91 


Hhere 


A^(cos  0-COS  0')^  f  4 


P  = 


A^icos  0-cos(d^  -t  (f)\^  +  4 


13.1.101 


If  the  annular  uir'oi!  is  at  nr-gle  of  attack,  then  to  these  induced  ve'cKMties  must  be 
added  those  induced  by  the  trailing* vortex  system  and  the  circulation  distribution  g^{Z).  The 
free-\ortex  system  from  the  duct,  when  it  i.s  at  an  angle  of  a'tack,  is  straight-line  tortices 
trailing  from  the  duct  to  infinity  and  parallel  to  the  2*axis.  Consequently,  they  do  not  induce 
any  axial  velocity  but  only  radial  and  tangential  velocities.  Those  induced  velocities  are 
given  in  Reference  15.  On  the  duct  itself,  the  contribution  to  the  axial-induced  velocity  due 
to  onglc  of  attack  arises  only  from  the  vortex  distribution  on  the  duct  and  is 


-  (.i;,.  ci,  Z)  ’  -  ^  /  “7“  r 

r)|.c/Z  '”j7,(Z)Jcos  ^ 


♦  3E[/c)-fC(k) 


13.1.111 


It  should  bo  noted  that  at  the  duct  leading  edge,  the  linearized  theory  gives  rise  to 
infinite  velocitie.s  unless  the  duct  scKi'tion  is  operating  at  its  ideal  angle  of  attack;  i.e.,  un¬ 
less  gjjfl)  •<  «o. 


4 


1 


1 


« 
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3.2  PRESSURE  DISTRIBUTION  ON  THE  ANNULAR  AIRFOIL 

Since  the  flow  field  is  nssumed  to  be  irrotational,  steady,  and  incompressible  and  body 
forces  are  neglected,  on  the  duct  itself  the  pressure  distribution  is  obtained  from  Bernoulli’s 
equation  as 


The  velocities  W^,  W^,  and  are  the  total  velocities  induced  by  the  various  singulari¬ 
ties  in  the  flow.  The  pressur  is  the  pressure  infinitely  far  ahead  of  the  duM  while 

Z)  is  the  local  pressure  on  the  duct  surface.  If  the  perturbation  velooUies  are  small, 
then  the  squared  terms  in  this  equation  can  be  neglected.  This  gives  the  linearized  Bernoulli 
equation,  which  on  the  duct  itself  is 


vr 

[3.2.21 

Since,  as  discussed  in  the  previous  section.  In  the  linearized  theory  the  axial  veloc¬ 
ity  induced  by  the  ring  vortices  has  a  singularity  at  the  duct  leading  edge,  so  does  the  pres¬ 
sure.  An  approximate  nonlinear  correction  to  the  theory  which  removes  this  singularity  can 
be  made  using  the  argument  given  in  Reference  15.  The  result  of  this  argument  is  to  multi¬ 
ply  the  axial  velocity  by  the  coefficient 
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13.2.3J 


As  normal  airfoil  shapes  have  infinite  slopes  at  their  leading  edges,  this  coefficient  will 
suppress  the  singularity  arising  in  the  velocity. 


3.3  FORCES  ON  THE  ANNULAR  AIRFOIL 


The  forces  on  the  annular  airfoil  follow  directly  from  the  Kutta-Joukowski  lawj* 
which  for  each  section  is 


F-pVr  tS.3.11 

The  velocity  V  is  the  velocity  perpendicular  to  the  direction  of  the  force  F  and  docs  not  in* 
dude  the  self-induced  velocity.  The  function  F  is  the  circulal.on  and  p  is  the  mass  density. 
To  obtain  the  total  force  in  any  one  direction,  the  force  on  each  section  is  integrated  around 
the  duct. 

The  total  lift  L  is  defined  in  some  srbitraty  direction  which  here  will  be  taken  to  be 
positive  radially  outward,  perpendicular  u>  the  axial  coordinate  a,  and  at  ■=  0.  For  the  axi* 
symmetric  case,  the  lift  on  each  section  is  radially  outward  but,  since  it  is  equal  at  each  an* 
gular  position,  *he  total  lift  is  zero.  By  the  sar.e  reasoning,  the  lift  in  an  arbitrary  axisym* 
metric  velocity  is  also  zero.  The  only  net  lift  occurs  when  the  duct  is  at  an  angle  of  attack, 
and  this  is  given  by^ 


fV  * - »  2w  fgAZ')dZ'=2n  f  g*(d)cos^  8d0  (3.3.2] 

/-/’‘■'t,  i  V  ^ 

The  induced  drag  also  follows  directly  from  F.quation  (3.3.1i.  As  for  the  lift,  the  in¬ 
duced  drag  <jf  the  duct  at  zero  incidence  is  also  zero.  However,  an  arbitrary  axisymmetric* 
radial  velocity  by  the  duct,  such  as  induced  by  a  propeller,  does  induce  a  drag.  The  induced 
drag  then  follows’  as 


For  the  duct  at  zero  incidence,  this  equation  reduces  to 


(3.3.41 
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where  the  velocity  is  the  arbitrary  axisymmrtfic-radial  velocity  by  the  duct.  Obviously, 
for  the  duct  at  zero  incidence  and  purely  axisymmotric  flow,  the  induced  drag  is  zero.  If  the 
duct  is  at  an  angle  of  incidence  in  otherwise  purely  axisymnietric  flow.  Equation  (.I.S.S]  be* 
comes,  after  some  manipulation. 


C„  =  —  [3.3.51 

which  is  the  result  obtained  by  Ribner. To  find  the  total  drag,  the  viscous  drag  must  be 
added  to  the  induced  drag. 

The  total  moment  on  the  axisymmotric  duct  at  zere  .ncidence  and,  also,  in  the  pres¬ 
ence  of  an  arbitrary  ax i symmetrical  flow  is  zero.  Consequently,  the  only  moment  arises  when 
the  duct  is  at  an  angle  of  attack.  This  moment  arises  from  two  sources,  one  being  from  the 
vertical  forces  (lift)  and  the  other  from  the  horizontal  forces  (drag).  This  differs  somewhat 
from  two-dimensional  airfoil  theory  where  only  the  lift  force  contributes.  The  moment  M  fol¬ 
lows  also  from  Equation  [3.3.1]  as 


M 
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~  Cu  +  Cl 
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where,  about  the  leading  edge 


t  .e.  0 


[3.3.71 


and 
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K(,k)-E{k) 


[3.3.81 


The  moment  from  the  drag  force  is  normally  small  compared  to  that  from  the  lift  force,  hut  it 
should  not  be  neglected. 
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To  obtain  the  moment  about  any  point  on  the  duct,  the  contribution  due  to  the  vertical 
forces  is  shifted  from  the  leading  edge  by  the  following  relation: 


f-i.- 


c^d-z) 


4.  COMPUTER  PROGRAM 

The  calculations  bast'd  on  the  theory  for  the  annular  airfoil  operating  at  a  given  angle 
of  attack  have  been  programmed  for  the  IBM-7090  high-speed  computer.  Input  consists  of  the 
section  camber  and  thicknc'^s  ordinates,  the  section  angle  of  attack,  and  the  chord-diameter 
ratio  of  the  annular  airfoil.  An  arbitrary  axisymmetric  velocity  can  also  be  included. 

The  standard  output  consists  of  the  lift,  drag,  and  moment  coefficients  as  well  the 
circulation  and  pressure  distribution  on  the  duct.  There  are  options  also  for  obtaining  the 
ideal  angle  of  attack  of  the  duct  section  and  the  velocity  field  inside  the  duct. 

It  takes  approximately  10  minutes  on  the  IBM-7090  high-speed  computer  to  calculate 
the  duct  forces,  pressure  distribution,  and' the  velocity  Held  inside  the  JucL  The  input- 
output  format  and  the  FORTRAN  listing  of  the  computer  program  are  discussed  in  th?  follow¬ 
ing  .sections. 


4.1  INPUT  FORMAT 

The  d'.<ct  identification  is  key  punched  on  the  fi  si  IBM  input  card  using  Format  12A6 
where  Columns  1  through  72  can  be  used.  The  remainder  of  the  input  data  is  p  inched  on  IBM 
cards  using  Format  F8.6  where  up  to  nine  field  point  p'.ra''  ''ters,  each  having  a  field  width 
of  eight  columns,  can  bo  punched  on  a  curd.  The.*-e  t  ticular  formats  are  de.scribed  in  Ref¬ 
erence  20. 

Seven  parameters  are  punched  on  the  s  co.id  input  card  with  the  first  four  quantities 
being  the  number  of  camber  ordinates,  the  nuinbe:  of  .section  thickness  ordinates,  the  chord- 
diameter  ratio  of  the  duct,  and  the  section  angle  of  attack  a  in  degrees,  respectively.  If 
the  camber  is  zero,  the  number  of  camber  ordinates  input  is  zero  and  if  the  ideal  angle  of 
attack  is  desired,  the  section  angle  input  is  zero.  The  number  of  camber  and  thickness  ordi¬ 
nates  which  cun  be  input  is  a  minimum  of  17  and  a  maximum  of  39. 

The  fifth  quantity  on  the  second  input  card  is  an  option  for  calculating  the  ideal  angle 
of  attack  of  the  duct  section.  If  the  ideal  angle  of  attack  is  not  to  be  calculated,  the  fifth 
input  quantity  is  zero;  if  this  angle  is  to  be  calculated,  ).0  is  used.  The  sixth  and  .seventh 
quantities  on  the  second  input  card  are  input  optinn.s  concerning  the  induced  radial  and  axial 
velocities,  respectively.  If  these  velocities  are  taken  as  zero,  the  input  value.s  are  zero;  if 
either  or  both  of  these  velocities  are  zero,  1.0  is  used. 


The  stations  along  the  section  chord  and  the  corresponding  camber  ordinates,  nondi* 
mensionaiized  on  the  section*chord  length,  arc  punched  on  the  first  and  second  sets  of  input 
data  cards,  respectively.  These  cards  arc  omitted,  however,  if  the  section  camber  is  zero. 

The  stations  along  the  section  chord  and  the  corresponding  half  thickness  ordinates, 
nondimensionalized  on  the  section  chord  length,  are  punched  on  the  third  and  fourth  sets  of 
input  data  cards,  respectively.  If  the  section  ideal  angle  of  attack  is  to  be  computed,  no 
other  input  data  is  required  for  the  computer  program. 

The  number  of  points  (a  minimum  of  17  and  a  maximum  of  39)  along  the  section  chord 
where  the  radial  velocity  is  pivon  input  should  be  punched  on  the  next  input  data  card. 
Following  this  card,  the  stations  along  the  section  ch..?d  and  the  corresponding  radial  '  cloc- 
ity,  nondimensionalized  by  the  ship  speed,  are  punched  on  the  fifth  and  sixth  sets  of  input 
data  cards,  respectively.  If  the  input  radial  velocity  is  zero,  th<-sfe  data  ca.xis  arc  omitted. 

The  number  of  points  (a  minimum  of  17  and  a  maximum  of  39)  along  the  section  chord 
where  the  axial  velocity  is  given  as  'nput  should  be  punched  on  the  next  input  data  card. 
Following  this  card,  the  stations  along  the  section  chord  and  the  corresponding  axial  velocit 
nondimensionalized  on  the  ship  speed,  are  punched  on  the  seventh  and  eighth  sets  of  input 
data  cards,  respectively.  If  the  input  axial  velocity  is  zero,  these  data  cards  arc  omitted. 

The  number  of  different  geometric  angles  of  attack  for  computing  th.>  duct  pressure 
distribution  (which  must  be  less  than  41)  is  punched  on  the  next  input  data  card.  Geometric 
angles  of  attack  in  degrees  and  the  corresponding  angular  positions  of  the  duct  section  in 
degree.^  at  which  the  calculations  are  to  be  made  are  punchr<d  on  the  ninth  and  tenth  sets  of 
input  data  cards,  respectively.  The  p:ogram  always  makes  calculations  for  °  0;  thus  it  is 
not  necessary  to  input  this  value. 

The  next  card  is  an  option  card  for  calculating  the  velocity  neld  inside  the  duct.  If 
this  velocity  is  not  to  be  calculated,  the  input  quantity  is  zero;  if  the  velocity  is  to  be  cal* 
culatcd,  1.0  is  used.  If  the  input  is  zero,  no  additional  r  ata  are  required  for  the  computer 
program. 

Two  input  quantities  should  be  punched  on  the  next  card  if  the  above  input  data  are 
1.0.  The  first  quantity  gives  the  number  of  different  radial  positions  (which  should  he  less 
than  or  equal  to  11).  The  duct  radii  and  the  stations  along  the  duct  chord  arc  punched  on 
the  eleventh  and  twelfth  sets  of  input  data  cards,  respectively.  An  example  showing  the  in¬ 
put  data  for  a  duct  where  the  radial  and  axial-induced  vclociiics  are  zero  and  the  pressure 
distribution  and  velocity  field  inside  the  duct  are  desired  is  shown  in  Appendix  A. 

Option  10,  which  allows  the  program  to  overload  the  operating  system.^”  must  also 
be  added  to  this  program.  The  binary  common  reassignment  card  must  have  a  minus  sign,  7, 
and  9  punched  in  Column  1,  a  plus  sign  punched  in  Column  2,  and  a  plus  sign  and  zero  puncl 
cd  in  Column  6. 
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4.2  OUTPUT  FORMAT 


Tho  input  is  printed  on  the  first  page  of  output  and  the  duct  forces,  including  the  lift, 
drag,  and  moment  coefficients,  arc  given  on  the  second  page.  The  third  page  of  output  gives 
the  duct  circulation  distribution,  nonlinear  corrections,  and  the  nonlinear  pressure  distribu¬ 
tion  inside  and  outside  the  duct  for  a  given  geometric  angle  of  attack.  It  should  be  noted 
that  the  pressure  distribution  includes  the  effect  of  adding  the  axial-induced  velocity  given 
as  input  to  that  induced  on  the  annular  airfoil.  If  more  than  one  geometric  angle  of  attack  is 
given  as  input,  the  parameters  presented  on  the  third  page  of  output  will  be  printed  on  a  new 
page  for  each  additionc!  geometric  angle  of  attac.*'..  The  axial  and  radial  velocities  along  the 
duct  chord  at  a  given  radius  inside  the  duct  are  printed  on  the  next  page  of  output  if  desired. 
If  more  than  one  duct  radius  is  given  as  input,  these  velocities  are  printed  or.  a  new  page 
for  each  additional  duct  radius  given  as  input. 

Note  that  if  only  the  section  ideal  angle  of  attack  is  desired,  this  computed  value, 
rather  than  the  duct  forces,  will  be  printed  on  the  second  page  of  output.  The  output  obtain¬ 
ed  from  the  input  date,  of  the  annular  airfoil  given  in  Section  4.1  is  also  shown  in  Appendix  A. 

4.3  FORTRAN  LISTING 

The  FORTRAN  listing  of  the  computer  program  is  given  in  Appendix  B.  In  addition  to 
the  subroutines  furnished  automatically  by  the  Bell  Monitor  System  on  thi  IBM-7090,  the  bi> 
nary  coding  for  the  BE-ELIP,  E2-AMGMHA.  B4-L.4CBRT.  DI-GLGAU2.  AMMATl,  AM-SIUF, 
VG-.AS+C,  ami  Kl-AQAKKl  subroutines  must  be  added  to  the  FORTRAN  listing  of  thi.® 
program. 


5.  EXPERIMENTAL  AND  THEORETICAL  RESULTS 
5.1  EXPERIMENTAL  RESULTS 

Two  ducts,  20  inches  in  diameter,  were  manufactured  and  tested  in  the  Subsonic  Wind 
Tunnel  of  the  Aerodynamics  Laboratory.  The  tests  were  run  at  a  Reynolds  number  of 
2.06  X  10^  Both  ducts  had  chord-diameter  rati  js  of  0.8  and  the  section  shapes  were  as 
shown  in  Figure  3.  The  cross  section  of  Duct  I  is  an  NACA  OCIO  thickness  distribution 
with  an  N.ACA  250  mean  line  of  maximum  camoer-chord  ratio  of  -0.0373.  Each  section  op¬ 
erates  at  an  angle  of  attack  of  6  degrees.  Such  an  annular  airfoil  typifies  the  shape  used  for 
Kort  nozzles.  Duct  II  has  a  cross  section  made  up  of  an  NAC.A  66  Mod.  thickness  distribu¬ 
tion  of  maximum  thickness-chord  ratio  of  0.10  with  an  NAC.A  a  *  0.8  mean  line  of  a  maximum 
camber-chord  ratio  of  0.04.  The  angle  of  attack  of  the  section  is  zero. 

Also  shown  in  Figure  3  is  a  duct  shape  te.sted  by  the  Bureau  Technique  Zborowski  in 
France.^*  This  duct  has  a  chord-diameter  ratio  of  0.96.  The  section  is  an  N.4CA  66-006 
with  no  camber,  and  tho  angle  of  attack  of  each  section  is  zero. 
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Force  data  in  coefficient  form,  i.e.,  lift  (C^),  drag  {Cp),  and  moment  (C^),  are  plotted 
in  Figure  4  for  Duct  I.  These  tests  were  carried  out  for  geometric  angles  of  attack  covering 
a  range  of  -  10  degrees.  Pressure  distributions  *  re  made  on  both  the  inside  and  outside  of 
the  duct  at  geometric  angles  of  attack  from  0  to  10  degrees  in  steps  of  2  degrees.  Flow  sepa¬ 
ration  occurred  on  this  duct  even  for  zero  angle  of  attack  thus  only  the  results  for  0  and  6  de¬ 
grees  arc  shown  here.  Figure  5  shows  (he  pressure  distribution  plotted  against  the  chordwisc 
location  of  the  pressure  tap  for  the  zero  angle  ol  attack  case.  The  results  for  the  b  degree 
angle  of  attack  are  plotted  in  Figures  6  and  7.  Figure  6  is  the  pressure  distribution  for  the 
duct  section  which  is  in  the  upward  position  as  shown  in  Figure  1;  i.e.,  ^  *  0  degrees.  Fig- 
are  7  shows  the  pressure-  distribution  at  th'i  lower  section;  i.e.,  6  »  180  degrees. 

The  force  data  for  Duct  II  are  plotted  in  Figure  S  and  again,  the  range  of  tests  was 
i  10  degret's.  Pressure  dtstributioi;s  for  this  duct  are  plotted  in  Figures  9  through  15.  Fig¬ 
ure  9  gives  the  pressure  distribution  of  the  inside  and  outside  of  the  duct  for  zero  angle  of 
attack.  Figures  10  and  1 1  give  the  pressure  distribution  for  4-degree  angle  of  attack  for  the 
position  ^  *  0  degrees  and  ^  =  180  degrees,  respectively.  Figures  12  and  13  are  the  pres¬ 
sure  distributions  for  8  degrees  and  »  0  degrees  and  <*>  >  l&G  degrees,  respectively;  Fig¬ 
ures  14  and  15  are  for  10  degrees  for  ^  a  0  degrees  and  ^  s  180  degrees,  respectively. 

The  force  data  for  the  BTZ  duct  are  shown  in  Figv*e  lt>.  •■r.d  the  pre.ssure  distributions 
are  given  in  Figure  17  through  19.  Figure  17  gives  the  pressure  distribution  for  the  zero  angle 
of  attack  case;  Figures  18  and  19  are  for  the  9-degree  angle  of  attack  with  ^  ^  0  degrees  and 
^  ■  180  degrees,  respectively. 

5.2  THEORETICAL  RESULTS 
5.2.1  Axisymmetric  Duct 

In  this  section,  results  of  theoretical  calculations  are  compared  to  the  test  results 
just  presented.  The  theoretical  calculations  were  made  on  an  IBM-7090  computet  by  the 
Applied  Mathematics  Laboratory  with  the  FORTR.AN  program  presented  in  Section  4. 

In  making  theoreticf.1  predictions  of  the  aerodynamic  characteristics  of  the  annular 
airfoil,  the  first  step  involves  the  calculation  of  the  section  circulation  distribution  as  given 
by  Equation.)  12.3.4)  anc  (2.3.15).  The  circulation  distributions  for  the  three  ducts  are  de¬ 
scribed  in  the  previous  section  are  shown  in  Figure  20  for  a  zero  angle  of  attack.  Since  the 
BTZ  duct  has  zero  cember,  the  circulation  di^-trtbutiun  is  solely  a  function  of  the  thickness 
tilciribution.  For  Duct  I  the  effect  of  local  angle  of  attack  is  ioininant  near  the  leading  edge, 
which  no  doubt  gives  an  indication  of  the  separation  which  occurred  on  the  duct.  For  Duct  II 
it  must  be  concluded  that  camber  is  the  dominant  effect  since  the  thickness  distribution  leads 
to  a  negative  circulation  distribution.  It  should  be  noted  that  these  circulation  distributions 
must  be  divided  by  y/l-z  in  order  to  obtain  the  actual  circulation. 


The  actual  circulation  distribution  is  singular  it  the  leading  edge  except  when  each 
section  is  operating  at  its  ideal  angle  of  attack.  This  ideal  angle  of  attack  is  obtained  from 
Equation  [2.4.3].  For  these  three  ducts,  this  equation  gives  the  ideal  angle  of  attack  as: 
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Desi^ 

degrees 

degrees 

Duct  I 

-1.67 

6 

Duct  II 

0.78 

0 

BTZ 

0.020 

0 

Tho  theoretical  pressure  distribution  on  the  ducts  is  given  by  Equation  [3.2.21.  Re* 
suits  of  evaluation  of  this  equation,  using  the  nonlinear  correction  given  by  Equation  [3.2.3] 
for  the  axisymmetric  case,  are  shown  as  .solid  lines  on  Figures  5,  9,  and  17  for  Ducts  1,  II, 
and  BTZ,  respectively.  The  prediction  is  poor  for  the  outside  of  Duct  1  where  separation  oc¬ 
curred  near  the  leading  edge.  Inside  the  duct,  however,  the  prediction  is  good.  For  both 
Duct  II  and  the  BTZ  duct,  the  agreement  between  theory  and  experiment  is  quite  good.  Also 
shown  in  Figure  9  is  the  pressure  distribution  for  Duct  II,  calculated  from  the  linearized  the¬ 
ory  a-nd  by  the  method  of  Chaplin.^^  This  last  method  is  a  nonlinear  thcoiy*;  therefore,  this 
plot  gives  the  comparison  between  the  linear  theory,  nonlinear  theory,  and  experimenL  It  is 
clear  that  the  liiiear  theory  does  not  give  quite  as  good  a  prediction  as  cither  the  nonlinear 
correction  or  the  nonlinear  theory.  However,  the  theoretical  pressure  distribution  by  any  of 
the  methods  is  quite  reasonable  except  when  separation  occurs  o.*;  the  duct. 

The  separation  occuiring  near  the  leading  edge  on  Duct  I  is  undoubtedly  a  laminar  sep¬ 
aration.  With  this  type  of  separation,  usually  called  leading  edge  bubble,  the  flow  reattaches 
itself  near  the  transition  region.  At  high  angles  of  attack,  turbulent  separation,  of  course, 
does  occur  and  the  anrular  foil  stalls.  These  various  flow  regimes  have  been  mapped  on  an 
annular  airfoil  by  Eichelbrenner.*^ 

In  addition  to  knowing  the  pressure  distribution  on  the  duct  (or  velocity  distribution), 
it  is  desira.ile  to  know  the  velocities  induced  anywhere  in  the  flow  field  by  the  presence  of 
the  duct.  These  velocities,  for  the  axisymmetric  duct  at  zero  angle  of  attack,  are  calculated 
by  the  comp-iter  program  from  Equations  [3.1.1]  to  [3.1.4].  Figure  21  shows  the  radial  varia¬ 
tion  of  the  velocity  at  the  1/2  chord  point  of  Duct  II  as  calculated  from  these  equations.  No 
experimental  re.sults  are  available  for  comparison. 
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5.2.2  Duct  at  on  AngU  of  Attock 


For  the  duct  at  an  angle  of  attack,  the  forces  and  the  effect  on  the  pressure  distribu* 
tion  are  independent  of  the  section  shape  except  for  the  moment  which  arises  from  the  hori' 
iontal  forces  The  forces  and  pressure  distribution  are,  of  course,  dependent  on  the 

chord>diameter  ratio  A. 

Figure  22  shews  the  theoretical  lift  coefneient  versus  the  chord^dirmeter  ratio.  Also 
plotted  on  this  curve  are  the  test  spots  from  Ducts  I,  II,  and  BTZ  along  with  test  results  by 
Fletcher. The  ducts  tested  by  Fletcher  had  Clark*Y  sections  with  thickness-chord  ratios 
of  0.117.  The  thi'ory  gives  a  surprisingly  good  prediction,  even  for  Duct  1  which  had  separa¬ 
tion.  The  theoretical  lift  is  also  plotted  on  Figures  4,  8,  and  16  as  a  solid  lino.  These  plots 
show  that  the  thecry  giie*:  good  prediction  of  the  lift  through  10  degree.-?. 

The  induced  drag  is  a  function  of  the  lift  sequared  as  shown  by  Equation  [3.3.5].  The 
theoretical  values,  calculated  from  this  equation,  are  shown  on  Figures  4,  B,  and  16.  Here 
the  drag  at  zero  angle  of  attack  was  taken  as  the  profile  drag  and  no  attempt  was  made  to  cal¬ 
culate  this  drag.*  It  can  be  seen  that  the  theoretical  drag  gives  a  good  prediction;  however, 
at  the  higher  angles  of  attack,  more  deviation  occurs  than  for  the  lift.  For  Ducts  1  and  II  the 
theory  underpredicts  the  drag  at  the  high  angles  of  attack,  which  no  doubt  is  due  to  separation. 
It  should  be  noted  that  the  much  higher  profile  drag  of  Duct  I  is  a  consequence  of  a  high  pres¬ 
sure  drag. 

The  theoretical  moment,  as  given  by  Equation  [3.3.6],  con.sists  of  two  parts,  one  due 
to  vertical  forces  Cu  and  the  other  due  to  horizontal  forces  C,,  .  The  moment  Cu  ,  as  given 

by  Equation  (3.3.7),  is  plotted  in  Figure  23.  Figures  4 ,  8,  and  16  show  the  total  theoretical 
moment  as  a  solid  line  compared  to  the  experimental  values.  .Also  shown  on  Figures  4  and  8 
as  a  dashed  line  is  the  moment  .  The  e^ect  of  the  horizontal  term  is  small  but  its  use 

does  give  a  better  prediction,  especially  for  Duct  II.  It  would  be  expected  that  the  moment 
of  Duct  1  would  deviate  somewhat  from  the  theoretical  value  since  the  leading  edge  consid¬ 
erably  alters  the  pressure  distribution. 

The  theoretical  pressure  distribution  contributed  by  the  angle  of  attack  as  calculated 
using  Equations  [3.1.11]  and  (.3.2.2)  is  shown  in  Figure  24  and  is  further  tabulated  in  Ta¬ 
bles  1  through  10.  It  can  be  seen  that  the  chord-diameter  ratio  has  a  considerable  effect  on 
the  pressure  distribution  and  that  the  shorter  the  chord,  the  greater  the  effect.  To  obtaiu  the 
total  pressure  distribution  on  a  duct  at  an  angle  of  attack,  these  values  are  added  to  those 
obtained  for  the  axisj-mmetric  case.  For  Duct  1,  the  pressure  distribution  is  shown  at  a 
6-degree  angle  of  attack  in  Figure  6  for  the  angular  position  ^  »  0  degrees,  and  in  Figure  7 


*The  ;aeU)o4  of  predtcUng  drag  of  asUyaimelric  bodiea  aa  ditcuased  by  GranviUc^^  ahould  give  a  reasonabl' 
eaUmate  of  the  profile  <feag. 


for  the  anguiiii*  position  <*>  «  180  degrees.  Because  of  the  flow  separation  occurring  on  the 
outside  of  the  duct,  the  prediction  is  poor  at  ^  >  0  degrees;  however,  for  ^  =  18U  degrees, 
the  predicted  pressure  distribution  is  quite  reasossbie.  It  ’  apparent  that  no  separation  is 
occurring  in  this  part  of  the  duct  at  the  6*degree  angle  of  attack. 

Theoretical  pressure  distrihufiens  are  shown  on  Figures  10  through  15  for  Duct  II  at 
various  geometric  angles  of  attack.  All  theoretical  pressure  distributions  look  reasonable 
up  through  8  degrees.  At  10  degrees,  however,  separation  occurs  on  the  inside  of  the  duct 
at  the  angular  position  °  180  degrees  (Figure  15).  At  8  degrees  there  is  some  deviatic‘> 
from  the  experimental  values  as  shown  in  Figure  13,  which  would  indicate  that  flow  is  start* 
ing  to  separate  at  this  angle.  The  predicted  pressure  distribution  on  the  outside  at  ^  »  0  de¬ 
grees  is  quite  reasonable  even  at  lO  degrees  (Figure  14). 

For  Use  BTZ  duct,  theoretical  pressure  distributions  are  plotted  with  the  experimental 
points  on  Figures  18  and  19.  The  geometric  angle  of  attack  is  9  degrees,  and  for  both  the 
^  c  O'degree  and  <f>  ~  iS0*degree  angular  positions,  the  comparison  between  theory  and  ex¬ 
periment  is  good.  Some  deviat  on  does  occur  near  the  trailing  edge  inside  the  duct  at 
=  0  dej-recs. 

CONCLUSION^) 

From  the  results  of  this  investigation  on  the  aerodynamic  characteristics  of  annular 
airfoils,  we  can  make  the  following  conclusions: 

1.  The  theory  gives  good  predictions  of  the  lift,  induced-drag,  and  moment  coefneients 
even  when  laminar  separation  occurs  at  the  leading  edge. 

2.  The  theory  gives  good  predictions  of  the  pressure  distribution  except  when  separation 
is  present. 

3.  The  moment  coefficient  from  the  drag  forces  is  normally  small  compared  to  that  from 
the  lift  force  but  should  not  be  neglected. 

4.  The  nonlinear  correction  to  the  pressure  distributions  is  normally  small  but  should  not 
be  neglected. 

5.  The  computer  program,  as  compared  to  other  methods,  is  efficient  and  quite  versatile. 

6.  The  computer  program  can  be  used  for  ducted  propeller  designs  with  the  restriction  of 
the  assumption  of  infinite  number  of  blades  since  an  arbitrary  axisv-mmetric  velocity  ran  be 
given  as  input. 

7.  The  computer  program  can  be  used  to  obtain  the  ideal  angle  of  attack  of  the  duct  sec¬ 
tion  which  is  important  for  design  purposes. 
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Figure  4  -  Lift,  DtAg,  and  Moment  Coefficients  for  Duct  I 


Figure  5  •  Pressure  Distribution  for  Duct  I  at  Zero  Angle  of  Attack  (a 


Figure  6  -  Pressure  Distribution  for  Duct  1»  a*  "■  6  Degrees  and  ♦  ■  0  Degrees 


on  for  Duct  I,  O-  «  6  Degrees  and  ^  >■  180  Degrees 


Pressure  Distribution  for  Duct  Il»  ct*  *  4  Degrees  and  e  ■  0  Degrees 


Figure  n  -  Pressure  Distribution  for  Duct  II,  >  4  Degrees  and  <»  «  180  Degrees 
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Figure  12  -  Preostire  Distribution  for  Duct.  11,  a^.  ■  8  Degrees  and  1»  «  0  Degrees 


lbiitif>n  tor  D<kI  II,  oi»  •  •  rw(T**«  anJ  «  ♦  t»«icrc«* 


flfm*  IJ  -  rr*«iur«  t^r  Ck.fi  t|,  •  19  «»4  «  •  iM  IK-tr#*. 
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Figure  20 


Radial  Variation  of  Velocity  at  lj|  Chord  Point  for  Duct  II 
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Moment  from  Vertical  Forces  as  Function  of  Duct  Angle  of  Attack  and  Chord-Diameter  Ratio 


Figure  24  -  Pressure  Distribution  as.i^Function  of  Duct  Angle  of  Attack  and  Chord-Diameter  Ratio 


Figure  24  -  Pressure  Distribution  as  .^Function  of  Duct  Angle  of  Attack  and  Chord-Diameter  Ratio 
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TOTAL  MOMENT  COEFFICIENf  ABOUT  L.  E.  (E03.3-6) 


ANO  X«5rATION<  along  CHORD  HHER6  Xa .S( |-COS( T ) ) 

STAR  aCIRCULATtON  OlSTHtSUTtON  tea2*J-tS| 

ICO  AND  NCIaNONLi»£AR  CORRECTION  OUTSIOC  AND  INStOC  DUCT  RESPECTIVELY  (E03«?-3> 

O  AND  PlaNONLINCAR  PRESSURE  DISTRIBUTION  OUTSIDE  AND  INSIDE  THE  DUCT  RESPECTIVELY  EOR  PAaO  DEGREES 
OP  and  pip  bmoNLINE^P  pressure  distribution  outside  and  INStUE  THE  DUCT  RESPECTIVELY  FOR  PAaINPUT 


t  ANO  XaSTATIONS  ALONG  CHORO  MHEflE  X**  %<  l-COSt  T )  ) 

CSTAR  rCtnCULATION  OISfRieUTION  (1:02.3-15) 

NCO  AMO  NCt«NCKLtNFAR  CORReCTION  OUTSIDE  ANO  INSIDE  DUCT  RESPECTIVELY  1E03.2-3) 
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ERASE  SUMt* SOM2,SUM3  1 

!00  1200  MbI.MPI  I 

SUMl»SUMl»SSM(Ml«FtM)  I 

SUM2=SUM2*SCM(M»1 ><SJT( tN.M)  f 

i  1200  SUM3=SUM3«SSM(M)«<.H(  IN.M)  j 

^  IFCATl  1316.1323.1316  I 

'/  1316  I  )>-2.*SUM24-SUM3  \ 

^  CO  TO  1300 

1 3.“  3  FK«  I  )«U2.*S1NF<  ALFA1/C0SF(AI.FA>«2.»SCH(  I  )  )*CTH(  I  I'FSTlU  1  )»»-2 

l.•SUM2♦SUM3) 1 

l.lOO  CONTINUE  I 

00  1301  K*t.MPl  I 

ERASE  FI .IN  j 

00  1201  l«tl.I3.I2  I 

tN3|N«  1 

1201  FlsFI«e{tN.K)*CTHC ll«a^«SM(ll«OX  1 

1301  Ftt(K)=FI  { 

00  lAOS  K«1 .MPl 

,  00  ISOO  Mst.MPl  I 

ERASE  OIM.FIH. IN 
j  OO  1400  tsll.!3*t2 

1N=1N41 

'  0tM=01H«B(  IN.K)*ST<  t  MSJTI  IN.MMSMI  I  )*0X  j 

1400  FlMsFtM>0(  1  N.K  lACBl  1N.P)*ST(  I  )«SM(  [  >*0.. 

I  OlMMIMlzOIM 

I  1500  FlNMIMlaFIM 

S4s2.«(SINFf ALFAt/C0SF(AUFAl4SCM( I) )»Ftl(K) 

;  ERASE  S1.S2.S3  ! 

00  1404  Mst.MPl 
SlsSt4SSM(M)*F(H)*Ft  I  CK) 

S2sS24SCM<M41 
1404  S3aS34SSM(M )«FtMM(M) 

IF(AT)1318. 1333.1318 
1318  A(KI=>2.*S24S3 
AIllsO.O 
GO  TO  1405 

1333  A(K)>S1~S24.S*S34S4 
1405  CONTINUE 

1F(AT)22C 9.2000. 2230  ! 

2200  MPlsMPl'l  j 

00  3333  Kal.MPI 
00  3333  3^1. MPl 
C(K.J)sC(K4 1.341) 

3333  A(K)sA(K4l) 

2000  CALL  MATINVIC. MPl. A. 1.0.10)  | 

IF(10-1>  1601.1600.1601 

1601  PRINT  CO  ! 

88  FORMAT! 14H  C  IS  SINGULAR)  ' 

CO  TO  1603 

1600  IF! AT) 1660. 1670. 1660  ! 

1670  Cl!  l)3FF!l)-.3ie309«A!l) 

IN^l 

19«11412  ■ 

00  17001al9. 13.12  ! 

IN«IN4l 

ERASE  SUM 

00  1701  3«2.MPI 

1701  SUM3SUM4.3I8309«STH!1)4A! J)«S3T!iN.J-l) 
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1700  Cl(  I  )«FF(  I  )-.318309*CTHf  1  )»A(  IMSUM 
CALL  SUB  AT 
29993  NH»36 

ERASE  NX 

DO  20001  I>tl«I3»12 
NX«NX^1 

20001  CX(NX)sCl(l> 

NNx2*NX-2 

KX«NX-2 

00  20002  KK«1«KX 

NXlsNX^KK 

NX2=NX-KK 

20002  CX(NXl)s«GXiNX2) 

CALL  GHHAS|NN.NH*GX*RY(200) ) 

AO-RV(2001 
00  20003  Jsl.NH 
JFIX*200-J«5 
AX( J)sRV( JFIX) 

20003  8X( JlsRYl JF |X>1 ) 

80=RY( 199) 

20004  FORMAT* 110.2E1S.6) 

ERASE  IN 

PRINT  103 

00  20006  t=ll«13*12 
IN3lN«l 
ERASE  Y1.Y2 
OO  2000S  Xsl.NH 
Y1=Y1*AX<K)»CJT< IN»K) 

20005  Y2sY2«eX(K) «SJT( IN,K) 

Y3=A0>YI>80^Y2 

IF( 1-1 )27000*2700l«27000 

27001  GG(I )«Y3/SINF< 1./S7«29b76) 

GO  TO  27002 

27000  GG< I >*Y3/STM( 1 ) 

27002  L-5  1-1 

XE=-< .5«( t.*CT( I ) ) )41 *0 

20006  CONTINUE 
ERASE  VI 

00  28001  i«i.iai«s 

28001  YlsYl»12<5664«Gl( 1 )•W< 1 )«CTHf 1 )*SM(I )«OX 
CDIsVl 

CT0=2.«M»C01/3«14t6 
28004  FORMATI 110.2CIS.6) 

26000  FORMAT (9F8. 6) 

26009  FORMATI 7E15. 6) 

PRINT  103 

ERASE  IN.J 

OO  1900  l«l U 13.12 

ERASE  XU 

INsIN^l 

OO  1901  '4«1.MP1 

1901  XU«XU»FLOAT'  < M ) •SSM(M) vcaTl I N.M ) 

1900  $U<I)*XU 
ERASE  IN 

00  20007  t8«11.13.t2 
tN«IN4l 

tF( 18-1 )200 09. 20008.20009 
20009  3SJ412 
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I 

i 

i 


( 

I 


2000a  TaFLOATFC J>/S7.29S78 

SLT»-CUBe*»TF<T  » 

OLT*CUBERTF( 3. lAIO-T ) 

FRa:CAUSS(SLr.ll|.r*NP.Or  ) 

AKaSQ9TF<4./(M«»2»CCTJ 10 )-COSF t DT»» 34T ) ) ••24A . ) ) 
IF(AK-1.)  20011.200t2.2.‘wat2 

20011  IF(AK>.99)  2001 A.200I 3.2001 3 
20014  CALL  ELLtPl Xt.AK.X3.X4.BF.BS) 

GO  TO  2001s 
20013  AK2sAK*42 

CALL  VSELLC AK2.bS»8F) 

GO  TO  2  0  OIS 

20012  ERASE  BS.BF 
2001s  ERASE  Y1.Y2 

00  2001 A  ICsI.NH 
2*teFLOATF<K»*IOT»»34T> 

Yl«Y14AX<K)«COSF(ZK) 

20016  Y2aY24BX(K>4SINF<2K) 

V3»A04Yt4804Y2 

PR»Y3*AK*<BF-BS)*C0SF< -S»COT«434T > ) •3.«OT»*2 
VG( lBI«.lS9tS4«H4PR 
00  1807  1x11.13*12 

AKxS0aTF(4./tH«424<CTl tB )~CT( 1 ) 1«4244. ) ) 
1F(AK>1*) 1802. 1810. 1810 
1802  IF; AK-.99)taoa. 1809. 1809 

1808  CALL  ELLIPI X1.AK.X3.X4.BF.BS> 

GO  TO  1807 

1809  AK2xAK*42 

CALL  YSELLIAKZ. BS.BF) 

GO  TO  1807 

1810  asxi.o 

1807  eVOl 1 )xSuf 1 )*AK«BS 
ERASE  GXO 
00  6666  1x1 1.13.12 
6666  GXOxGXO^BVOll) 

X  ASM! l)»0X/3. 1416 

00  6662  Mxl .MPl 
ERASE  GXX.INN 
00  6661  1x1 1.13.12 
1NN«1NN41 

6661  GXXsGXX42.«BV0f 1 )«SMI 1 )4CJTf I NN.M > aOX/a. 1 4 1 6 

6662  GX(M)sCXX 

ERASE  SUM 10. SUMO* SUMS 
00  11905  Mxl. MPl 
SUMt0xSUM104GXfM)*SJT( IN.M) 

SUM9XSUM9-FL0ATF ( M ) •CX ( M ) *2 • 

1 1 905  SUMSxSUMS'i'Z  .  •  (  - 1  •  x^M  )  XFLOA TF  (  M  )  XGX  t  M  ) 

IF(lB-l)  11912*11908.11912 

11912  1F<I8-181)  11907*11906.11907 

11907  «0x<-2./STf IB>)*(GX04SUM10) 

GO  TO  20007 

11908  VQxSUMO 

GO  TO  20007 

11906  W0XSUM8 
20007  VOl lB)x«0 

ERASE  IN 

20000  00  20020  18x11.13*12 
INxlNxl 
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tF( IB-I >20021 .20020.20021 
I  HV6sVG(  1B)'»^T2(  IB) 

C1STH3.5«GI ( IBl/STHC IBl 

•G0®MV6-G1STM 
•GI=HV6«G1STH 
'  woa=  VOCIB) 

IIGG=VG(  IB) 

PO'=2.«(  WGO^WQO) 

P|^2.*OlGI»MOO> 

POL* IB)=PO 
PILI IB)=fM 
PON=S( IN)*P0 
P|N=y<  IN) •P  I 
X( IN)=.5*P0 
XX( IN)=.5*P! 

}  CONTINUE 

]  format ( I  10. 7e IS. 6) 

GO  TO  1603 
}  ERASE  IN 
A( I )=0.0 

00  1669  1=11.13.12 
iNsiN'fl 
ERASE  SUM 
DO  1661  J=2.MP1 

1  SUM=SOM»A<J )*SJT( IN. J-1) 

GUI  )  =  FF(  1  )  4.  tS91S««SUM 
CGI I)=S1( 1) 

IX=I-I 
9  CONTINUE 
ERASE  T.IN 
00  1234  I=I 1.13.12 
IN=IN4l 

4  T=T  481  IN,  1)  *0  I  1  )4ST(  I  )4SM(  I  )«0X 

AnTsATANFI . 079577*T-.54SUM1-SCM( 1 )) •57.29578 
ERASE  AT 

3  FORMAT (55MI lOEAL  ANCLE  OF  ATTACK  IN  DEGREES 
1F9.6) 

PRINT  1753. ANT 

O'J  20024  1=1  1.13.12 

4  Gl  I  n^GlI  1  )  •STmI  I  ) 

IF(AT)  1667.1668.1667 

7  CALL  SUB  AT 
CO  TP  29993 

2  FORMATI 1 10.5E15.6) 

1  F0RMAT(8F9.6) 

0  FORMATISIA.AFIO. 6/1814) 

2  FORMAT(eF10.6) 

1  FORMAT! 8E 15.6) 

3  CALL  SUBmR 

8  CALL  END  JOB 

eAp 
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2901 


21993 

21990 

21998 


2601 


2599 

2602 


2598 

3000 

21005 


FOR 

SULROUTINE  SUB  HA 

DIMENSION  SSM( lei ) .SCM( 1 81 > * CT< 18 1 ) . 5T( 1 8 1 ) fCTHl 1 81 ) .STH( 101 ) »SJT( 
140  1 40)  .CJT( 40.40) . 01 4 0.20 ) . SM ( 1 82 ) * APM( 20.20) .F < 20 ) . CB( 40 . 20 ) .8( 40 
2.20).C(20.20) .FF( 1 6 1 ) . A( 20 > . Gl ( 1 8 1 ) . W ( 1 8 1 ) • SU C 1 Q 1 ) « VG ( 200 ) . V0(  200 ) 
3.  bVO(  181  ). 

4FII(20).  OIMM120).  FtMM(20).AX(  100)  .0X(  100).HY(200).GX(200) 

COMMON  SSM, SCM.CT.Sr .CTH.STH.SJT.C jr .G. SM . APM .F .CP. O. C .FF . A . Cl .W.S 
lU.VG.VO.HVO  .'r  t  I  .OIMM.FIMM.AX.HX.HY.GX.H.  I  1  .  1  2  .  I  3  «  OX  .  X  1  .  X  3  .  X4  .  NP  .  NH 
DIMENSION  Y ( 161 ) . a( 181 ). TITLE! 12)*X(101).XX(lPt).XC<40.2).XT(^0.2) 
l.XSC(40.2) 

COMMON  Y.S. T I TLE.X.XX. XC.XT .XSC 
D I  MENS  I  ON  AXl ( 4U ) .HXl! 40) .AX2(40 ) .UX2(40) 

COMMON  BF.8S.AK.A01.A02.ALFA.MP1 

DIMENSION  oG(  181  )  .GP(  lEl  )  •  SMZ 1  20  )  .  X  I  (  20  )  .  T  I  (  1  8  1 )  .  T2 1  1 8 1  )  .  ZU  20  )  .  Z 
12!Z0) . WR2(20.40) .WT2(20.40) 
common  GG.GP.SMZ.XI.T1.Z1.Z2.T2.MR2.wT2  .WK  .AT 
DIMENSION  WR4(  20.4  0).WT4(20.AO).WRt>(  20) 

COMMON  MR4.Mr4.MR6.-C0I 

DIMENSION  POL! 181).PIL(181) .POPL! 181 ) .PIRLI IB  1) .POO! 181 ) .PI  I ! 181 ) 
COMMON  POL.PIL.PORL.PlRL.POO.pl l.COl .CTD.ALFAR .CL.CO.CM.CMl .CM2.PH 
1  1 

ERASE  IN.T1.T2 
DO  2901  i>i.iei.s 
lN=tN4i 

XX! 1N)3.S*!  1 .-CT!  I  )  ) 

X! IN)=FLOATF! I-) ) 

PRINT  103 
ERASE  XC.XT 
ERASE  AX. ex 
HEAD  21993. TITLE 
PRINT  2 1990. TITLE 

READ  219</e.CC.  TT.H.  ALFA.AT.WK.MW 

PRINT  1752.CC.TT.H.ALFA.AT.WK.WW 

NC=CC 

NT-TT 

M=NC/2 

FORMAT! 12A6  ) 

FORMATI IHJ.  12A6) 

F0RMAT!9F8.6) 

ERASE  XSC.SCM 
1F!CC)2601. 2602.2601 
READ  21998. !X! I ) . 1=1 .NC) 

READ  21998.  !Y!  I). 1  =  1. NO 
PRINT  1752. !X< 1 ). I=1.NC) 
print  1752.  !  Y!  I  )  .  1  =  1  .NO 
DO  2599  1=1.37 
S1=XX( 1 ) 

CALL  OISCOT !Sl .S I . X. Y. Y.- 120.NC . 0. S2 » 

Tl! 1 )=S2 

READ  21998. (X! 1 ) .1=1 .NT) 

READ  21998. ! Yl I ) « I=1.NT) 

PRINT  1752. !XI 1 )« 1=1.NT) 

PRINT  1752. !Y!  1  ).1  =  1.NT) 

DO  2598  1=1.37 
SlsXX! I ) 

CALL  DISCOT !S1 .SI .X.Y.Y.-tZO.NT.O.SZ) 

T2! I )=S2 
format !3E 13. 4} 

FORMAT!  I  10.  f'ElS.6) 


T2 


000  FORMATC  Ia.«FfS.^  ) 

00  218B8  IC=i.2 
00  2597  1=1,37 
597  Y(  I  )  =  T1 ( n 
Y<  1  )  =  ,0 
■  Y<37>=.0 
QOl  F0RMAT(9E13.4) 

498  FORMAT(3E15«6) 

00  21001  KK  =  1,.IS 

NX1=374KK 

NX2=37-KK 

001  YlNXl )=-Y<NX2) 

IF(1C*2)  2408,2409.2408 

408  IC=2 

CO  TO  2410 

409  1C=1 

410  00  21002  IS=1,  IG 
1F(T1(51>  21977.21978,21977 

978  ERASE  RY 

CO  TO  21976 

977  CALL  CRHASl  72,36. Y.RY( 200) } 

976  AO=RY(200) 

XC(l,tS)=A0 
00  21997  J=1.36 
JF 1X=200-J45 
AX< J)=HY( JF IX) 

XC(J41, 1S)=AX(  J) 

BX(3)=RY( JFtX-t) 

997  XT(3.1C)=ex( J) 

BOsRYl 199) 

ERASE 

996  FORMAT! I  10. 2E1S,6) 

ERASE  IN 

00  29000  1= 1, 181.5 
IN=!N41 

ERASE  Y1.Y2.SS 
L«I-1 

00  21995  K= 1.36 
Yt=Yt4AX<X)«CJT< IN.k) 

Y2=Y24ex(K) vSjTI IN.K) 

1F(  15-2)24000.23000.24000 
1000  erase  S 

CO  TO  21984 

>000  tF(  I-t 121983. 21982. 21983 

1983  1F(  !'*181)2198t. 21980, 21981 
1982  T0Ss2./SINF( 1./S7. 29578) 

GO  TO  21995 

1980  70S=2./S1NF ( 1 79./S7.29S78) 

CO  TO  21995 

1981  T0S=2./ST( 1 ) 

1995  SS=SS-FLOATF{K  )4BX(F)*C3T< IN.K)4T05 
YC  1N)=-SS 
XSCl IN, IC)=-SS 
SI  I )=-SS 

1984  Y3=A04Y14B04Y2 
>000  CONTINUE 

00  21006  K>;  =  1.3S 
NX1=374KK 
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NX2*3/>KK 

21006  YlNXl )»4Y(NX2) 

Y<3r»»Y<36) 

21002  V(36)3YC36) 

00  2596  I«l«37 
2596  T1(1)>T2(1) 

21688  CONVINUE 

21994  FORMAT! 110. 2E 16*6* 1 10* 3E 15.61 

ALFAsS I NF(AUFA/S7.29578>/C0SF( ALFA/57. 29578) 
29001  OO  22001  t=1.^7 

SC  I  )st  ./SQRTFC  1.4lXSC<  1.  1  )4ALFA4XSCC  I  .-ri  )«  '2) 
«'200t  Y(  I  )  =  1./SQ9TF(  t  . ♦  ( XSC  C  1 .  1  )  4 ALFA-XSC  (  1  *2f  )*42) 
00  22889  J=t.36 
SCM( J)s-xC( 3.2) 

22889  SSMC J)sXT< J.2) 

SC  1 )3.0 

YC  n=.o 

SC37)s.O 

YC37)a.O 

ERASE  IN.W 

00  22222  Isl.lal.S 

tN«IN^l 

X( lN)aFLOATFCl-l) 

22222  XXC IN)=.S«C l.-CTC 1 ) ) 

ERASE  IN 
|FC«K)96. 97.96 
97  ERASE  ■ 

CO  TO  190 
96  READ  21998* ANR 
NRsANR 

REAO  21998. CTtC 1 ).I«t.NR) 

REAO  21998. CT2C1). 1-1. NR) 

PRINT  17S2.CTlCl).|st.NR) 

PRINT  1752. CT2C I ). l^l.NR) 

ERASE  IN 

00  2699  1«1. 181.5 

INa|N4t 

SlsXXC IN) 

CALL  O1SCOTCS1.S1.T1.T2.T2.-120.NR.0.S2) 

2699  DC  I )«S2 
190  IFCWV) 196.197.196 
197  ERASE  T2 
CO  TO  95 

196  NEAO  21990. ANA 
NAxANA 

REAO  21996.CT1CI).1>I.NA) 

REAO  21998. CT2C 1 ).1«1.NA) 

PRINT  l7S2*CTlCl).t«l.NA) 

PRINT  17S2»CT2Cl).l:xt.NA> 

ERASE  IN 

00  2698  Isl.ISl.S 

tN*-IN4l 

SIsXXC IN) 

CALL  OlSCOr CS^.Sl.Tl .T2.T2.'120.NA.0.S2) 

2690  TIC  I  )  =  S2 
1752  F0RMATC9F10.6} 

103  FORMAT C IMI) 

95  RETURN 


01 MENS > OH  SSM( 181 l.SCMC 181 l.CTl leil.src 181 l.CTHl 181 ) .STH( 181 l.SJTl 
140«40) .CJT( 40.40 )tG( 40 tZOl .SM( 1 82 ) . APM( 20*20 ) .P < 20 > .CB ( 40 . 20 ) . B( 40 
2.20>.C(20.20).FF( 181 ) •A120>.G1( 181) .M( 181 ).SU< 181 ) . VG( 200 ) . VQ( 200 ) 
3.  BV0tl81). 

4F)1(20>.  01MH(20).  F 1MM( 20 ) . AXf 1 00 ) .BX< 1 00 ) .RYl 200 ) .GX ( 200 ) 

COMMON  SSM.SCM.CT.ST.CTH.STH.SJT.Cjr.G«SM.APM.F.CB.9.C,FF.A.Gt .W.S 
1U.VG.V0.BV0.fi  1  .0IMM.FIMM*A;.*BX.RY.GX.H.1  1.12.1  3cOX.X;.X'^.X«.NP.NH 
DIMENSION  Y( 181).S( 181 )• TITLE I  12) .X( 18 • « . XX| I  81 ) .XCC 40. 2 } .XT (40. 2 ) 
l.ASC(40*2) 

COMMON  Y.S.TITLE.X.XX.XC.XT.XSC 
DIMENSION  AXl(40).aXl(40).AX2 140) *8X2(40) 

COMMON  BF.BS.AK.A01.A02. ALFA. MPl 

DIMENSION  GG( 181 ).GP( 181 ) * SMZ( 20 ) .X 1 ( 20 ) . T1 ( 181 ) .T2( 18 1 ) .Z 1 ( 20 ) . Z 
i2(20).»R2(20.40).WT2(20.40) 

COMMON  GG.GP*SMZ.XltTl*Zl.Z2.T2.«R2.«T2  .«K  .AT 
DIMENSION  «R4(20.40).vr4(20.40).«R6(20) 

COMMON  WRA.WTA.WRA.rni 

DIMENSION  POLl 181 ) *P1L( 1«1 ) .PORL( 181 ) .PIRLf 181) .P00( 1 8 1 ) *P1 1 ( 1 81 ) 
COMMON  POL.PIL.PORL.PIRL.POO.PII.COI.CTO.ALFAR.CL.CD.CM.CMI .CM2. PH 
IS 

DIMENSION  SCO( 181).CXT( 131) .P( 181 ) 

DIMENSION  R( 181) 

ERASE  PI t*POO.ALFAR*PHl*CL*CO.CM*CMl.CM2 

ALFAlsALFA 

NM«10 

NPP«$ 

KlslO 

31*10 

NH*36 

ALFAR*1. 

PHls.O 

SS6  ALFAsALFAR 
)94  FCPHAT  (8E1S.6) 

198  FORMAT (OF?. 6) 

X1*0. 

ERASE  IN 

DO  11  1*11. 13.12 

IN«fN4l 

93  DO  10  K«1.NM 
M*K-1 
NP*NPP«M 

PP'.cGAiJS3(0.0.3.1A16.NP.CT) 

CMC*COSF(3T)  -CT(1) 

.F(CMC)99.2.99 
99  AK*SQRTF(4./(M 

1  *»24C  tr4»244.  )  ) 

IF(  A.'t  -1.  )11 

til  lF(AK-.V9)t08.109* 109 
108  CALL  E;  LIP(  XI. AK.X3. X4.BF.es) 

^.OTO  107 
<09  AK2«AK4*Z 

CALL  Y£€LL(AK2*BS*BF) 

GOTO  107 

2  2*0. 

GOTO  80 


I 


107  Z»<  (  (2.>AK««2)»*2»BS-4**(  l.>AX*«2)*»2*BF)/l'CMC«AK*«3)-I</CMC) 
60  IF(M)6.6.5 

5  2«2«C0SF<FLCATFIM)*0T) 

6  PRs2 
IF  (M)7*/.8 

7  8(  tN«KI=-l.S70a«H't-.3tS309»PR 
GO  TO  10 

8  8(  IN.K)  =  .63o(>ie*PR 

10  CONTINUE 

1 1  NpxNPP 
20  FORMAT (  1814) 

40  FORMAT! l4/( )CF7. 3) ) 

0X=FL0ATF<I2  )/171. 68734 
00  SOI  1«/.13.2 
SM(I.*ts4.0 
501  SM<l41)s2.0 
SM( 1  )«1.0 
SMt 13)=1.0 

8000  ERASE  C 1 1 .F 1 .C 12 • IN 

CO  1100  1  s  II  •  13.12 
2NZIN41 

Cl 1=C1 63061 e«B< IN. 1 )«CTHt I )*42*sm( Il^DX 
Fl=Fl-2.*ALFA*e<  IN.  1  )«CTH(  I  )««2*SM(  1  )  •D>:/'57. 29578 
1100  C123C124.318309*B( IN. 1 )«ST( I )«*2*SM! 1 )«OX 
C< 1.1)=C1I 
C! 1«2)=C12 
F(1)*F1 
00  100  3^3. J1 
ERASE  XCl.lN 
OO  1200  1=11.  13.12 
tN3lN4| 

1200  XCl3XCl«-.318309*ST(  1  ) *5  JT(  I  N»  J- 1  )  •SM t  I  )*OX*OC  IN.  1  ) 

100  C5  l..t)=xCl 

00  200  XxS.Kl 
ERASEXCKl  «XF 
ERASE  IN 

00  1300  1=11.  13  .12 

1N=IN4| 

XCKI  =  XCK1-. 63661 8«8( IN,XI*CTH< I )4»2  4SM(  I  )*DX 
1300  XF=XF+2.»AL  A4BC IN.K » •CTH( I > ••Z«SMC I ) *0X/57.29578 
F(K)xXF 

200  C(K, 1 )=XCK1 

OO  ISOO  K  X  2  .  K1 
00  ISOO  J  X  2  .  31 
*-RASE  XC.tN 

UO  1600  I  =  11  .  13.12 
IN=1M41 

1600  XC=XC-. 31830948! ir  ' 5 *S3T! 1 N. 3-1 ) ^ST ! I ) •SM ! I ) 40X 
1500  C!K.3)=XC 

OO  700  Kxl.Kl 
00  700  3x1.31 
IF  !K-3)702.701 .702 

701  C!r..3)=l.-C!K,3) 

GO  TO  700 

702  CIK.3)=-C!K.3) 

700  CONTINUE 

00  799  Kxt.Kf 
799  A!K)xF|K) 
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FORMAT  (///////eel 5. 6) 

I  CALL  MATlNV(C.Kt.A. I.X.IO) 

IF  ( tO~t )90l*900.901 
PRINT  ee 

.  FORMAT  (14H  C  IS  SINGULAR) 

GO  TO  603 

I  ERASE  SUMA»SG«SG1«NX* IN 
00  1000  I>1. 13.12 
INcIN*-) 

SUMA  >  0.0 
00  1001  J=2.Jt 
SUMAsSUMA^AI J)<SJT(1N. J-1 ) 

1F(  1-1)  1003. 1002.1003 
!  R( 1 )s( 2.«AL^A/S7.29578)-.636618«A( 1 ) 

GOTO  1004 

I  R( 1 )=2.«(ALFA  /S7.29S78)«C1H( 1 )  -.63661 8«STH( I ) • ( A( 1 ) •CTH( 1 > 

1  /STH(l)  4SUMA) 

»  SG=SC4R(1)*CTH(I)«SM(1  ) 

NX=NX4 1 
GX(NX)=R(t) 

I  SG1»SG14R(1  )«STH(  1)**2«CTH(  I  )*S:«(  I  ) 

CL=6.2e32*SC«OX 
CM=-6.2a32*SCl •OX 
C03.079S7r*H«CL«*2 
i  ERASE  CM1.CM2 
I  ERASE  FI. IE 
00  2497  1*1 1.13. 12 
CXT(1)=C)(I )*CTH(I) 

'  Fl3F14Gl(I)«CTH(l)«SM(l)«OX 
00  2200  in  1.13*12 
ERASE  XSCO 
00  2100  !P=1 1.13.12 

AK*S0RTF(4./<H««2«<CT( I )-CT( IP)  )«*244.  )  ) 

IF' AX- 1 .IS) 11.2100.3111 
I  IFfAK-. 99)3108. 3109. 3109 
-  )  CALL  ELLIPI  X1.AK.X3.X4.BF.BS) 

GO  TO  3107 

>  AK23AK*«2 

call  YSELL< AK2.BS.BF) 

f  XSCO=XSCOA(R(IP)«CTM( IP ) •H» C CT ( I )-CT( IP) )/AK) • (BF-BS) •SMI  IP) AOX 
)  CONTINUE 
)  SCOdIxXSCP 
ERASE  CHI 

00  2300  I«{ 1.13.12 

>  CMlsCMdCXT  (  1  )«SCO(  1  )^SM(  I  )«OX 
CMls'».25^CL^Fl«CMl 
CM2xCM4CMl 

NNs2*NX-2 

KX>NX-2 

00  2301  KKxl.KX 
NXlsNXAKK 
NX2XNX-KK 
1  GXINXI )«-GX(NX2) 

CALL  GMHAS(NN.NH.GX.RY(200) ) 

A03RY(200) 

00  5000  3x1. NH 

3FIXx200-3^S 

AX(J)xRV(3FIX) 
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toeo  •«( ji««TC jF iMoi » 

•o«itvc 

600*  roAMAf ( I 

OO  50«J  l>l 
CMASC  VI .Tt 
L-l-l 

OO  0003  K■I•NN 

TJC-FLOATr  (L  •«  t/%r*2*sro 

VlaVt»AXtKI •COS^f tK> 

5003  v»«v7*ox<Kt«siNrti«> 

T3«A0*YI«II0*V> 

5093  CONVtNue 
eilAS€  IN 

OO  7500  lal I•(3•t7 

Ca^LOATFI S- t )/57«7957A 
SLTa-CUOeOfFIC I 
UCTsCUH^ATF ( 3*  tAIA-C« 

OOaCAUSS  (  SL  f  »«A.  f  •  NPi*  •  OT I 

AKaS0nTF<4./<M««7*<CT( I »-<OSF|OT»»3*C I  I AAt^A. 01 » 

IFIAK'I.O)  7502»75t0*75t0 

7502  2FIAK>««9t  750O* 7509* 7509 
7505  CACU  CLLtFI K|.AR*a3««A.OF«9S> 

CO  TO  7507 

7509  AR2«AK«a2 

CALL  VSCLLI AK2«as*ori 
CO  VO  7507 

7510  ERASE  eF«B3 
7507  EMAuE  VI .VZ 

OO  7503  R«I*N*« 

ZRsFLOATFl*.  |•COT••3•E} 

VlariAAXlR) •COSFIZK) 

7503  V2aV24ei(Rt •SINT  1 ZKl 
V3>AOAVI«00«V2 

0««<0T««2»V3/AR1«<  <A*-A«*»2l»BS-l A.-3«AAK**2I *5?  »*COS7 l*5*10T**l*t 
1)1  •10T«A3*E»*«2 
Filial- 

X  3.*H/3«lAI5laP9 

1F( l-t 15701 •0700*0701 

5700  PIK  ll-*0 
P001l}««0 
CO  TO  7500 

5701  POOl  I  lalPIl »-R< I I/STmI 1 1 laCOSOEIPMI I 
PIICI|a{P<l)aR<l  l/STHt  t  )  laCOS^tPHt  I 

7500  CONTINUE 
2499  FORMATI I 10« 7E15«5I 
50  F0RMAT19r«.6l 
7555  ALPAaALPAI 
RETURN 

503  CALL  ENOJOe 
CNO 
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rom 

sv«»ouTiNC  su(i«n 

SSHI  ISI) •%€•*(  I •!)•€?(  l8t»*STI  tAl  I  •CThI  1 8t  )  •STM(  t at  )»SJT( 
t*0*«OI«CJt(«0*«0>«Ct«0.70l*SMl ta2t.*i>M(70.20).f t 20) •Cal«0t20).Of «0 

2.201  •ci.T0*2O).rr(  tat >.*(20t.&i( laii.ai lai i.soi lai  »•  vc(200) .vo(20u.' 
2.  avoitait. 

«rtii20i.  oi««(20i.  r tool.axt iu0>.«v(200t.cx(200) 

COMMON  SSM.SCM.CT.Sr.CtM.STH.SjT.Cjr.C.SM.APM.f .ca.B.C.FF. A.Ct .W.S 
tU.VC.VO.avQ.f I I.OtMM.A IMH.AX.ax.RV.OX.M. t  t . t2.l3.0X.Xt.X3.XA.NP.N*4 

otMCNSioM  VI  tail. St  tat  t.t itcei t2i.xt lat  t.xxi lat t .xc(40.2>.xTtAC^2> 

l•XSCI«0•2} 

COMMON  V.S. riTLC.X.XK.XC.XT.XSC 
OlMCNStON  AXtiAOl  .(tAt  (AO  t  .4X21401.0X21^0) 

COMk.ON  aA.aS.AK.AOI.A02.A(.FA.MPl 

oiMCNSioN  cci tai ).cp( iOi ).SMC(20).xt(20).ri( tat ).T2( lai ) .7t (20).z 
t2120).««2(20.40).Mr2(20.40l 
COMMON  CC.GP.SM/.Xl.r 1.2 . .22.12. •M2.Mf2  .MX  .AT 
OIMCNSION  •R4(20.40).Mr4(;  0.40).M<><>(20I 

COMMON  MMO.MTA.MRa.COI 

OIMCNSION  Roci lat I.RIL. lat i.PORLi tai i.PiRLi tat  t.poot  tat  >.Pii ( 161 1 

COMMON  POt.Ptl..PORt..P|Mt..P00.Ptl.CUl.C1O.ALFAR.CL.CO.CMtCMt.CM2.PH 

II 

oo  iiaol  iBiuti.is 

CH4SC  ORC.aMO 
oo  Itaor  |Pall.l3.t2 

AK«soMrr(4./tM««24irT( t i-cti ipi >44244.) > 
ir(AK>t. )iiao2. I lato.t laio 
11002  iriAKo.ooti taoa.ttaoo.iiaoo 
Itaoa  CACL  Cl  L|P(  Xt.A«.X3.X4«eC.»S) 

CO  TO  tiait 

Itaoo  AK2«AK442 

CALL  VSCLtI AX2.BS.e2t 

CO  TO  t tail 
ttaio  CMASe  BRC.ORQ 

nait  aM0«ORO<'.3tO30«4H4SU<  (PIaIBP-BSMSMI  |P)40X 

I  taOT  eRC«aRC4.029S7  24Ct ( tP)«AX4(H*424(CT( t )>CT( IP) )  ••24(8F>US)>2.4ttS) • 
tCTHt |P)4SM( tPtVOX 
VOI t laBRO 

tiaol  vc(i)«eRC 

ERASE  tN 

00  12006  lall.lS.IE 

INatNAl 

L»l'l 

XEao(.S4| l.ACTi I  I ) >41.0 
aOOaVOI I I-C.VSUI I ) 

VOI«VO(  1  )«'2.«SU(  1  ) 

MOaVCI I l♦«0O 
Mlavct  t)4*Ot 

POHa«; l-t .4X1  IN) )442«t .4X0442) 

PORN«S( IN)4P0R 

PIRa-U-t.4XX(  INt  )••^-^.4»^44^^ 
l>|RN«V(  1N)4PIR 
PORLI 1 )aPOR 
PIRLI 1 )aP|R 
12006  CONTINUE 

1753  20HMAT(55H  CHORD  DIAMETER  RATIO 

1C9.4///55M  LIFT  COEFFICIENT  PF-  DECREE  ICOT.3-2) 

)F9.4///55H  induced  DRAG  COcFFlCIENT  (£03.3-5) 

1F9.4///55M  MOMENT  CCEFFICIENT  FROM  VERTICAL  FORCES  tEQ3.3-7) 


IF9.4///S5H  MOMENT  COEFFICIENT  FROM  HORIZONTAL  FORCES  (C03.3-a) 
1F9.4///SSH  TOTAL  MOMENT  COEFFICIENT  ABOUT  L*  E.  (E03.3-'6) 
1F9«4///SSH  iNOUCeo  DRAG  OUE  TO  RADIAL  INDUCED  VELOCITY  IE03.3-4) 
IF9.4) 

PRINT  I753«H«CL«CO»CM»CMl,CM2.COI 

READ  26000. AXA 

NOA«Ar-;A 

READ  2C025. (AX2(lAI.IA=t.N0A) 

READ  2002S. f AXlI lA). |A=l.NOA) 

DO  8809  lAsl.NOA 

trSS  format  (S5;il  GEOMETRIC  ANGLE  OF  ATTACK  IN  DEGREES 

1E13.4  ZSSH  POSITION  OF  SECTION  AT  ANGLE  fPA)  IN  nCGREES 
1E13-4/I 

PRINT  t7SS.AX2(  lAI.AXM  lA) 

1756  FORMAT! ItSH  T  X  GSTAR  NCO  NC 

11  PO  PI  POP  PIP  ) 

PRINT  1756 

AXl f |A)sCOSOF( AXlf IA>  > 

ERASE  IN 

DO  7777  lsll.13.12 
IN»1N41 

POT>(POL< 1) APOOl I }«AX2C lA) >*SI IN) 

PlT«(PIL(l}4Plt< t)«AX2CtA))*VIlN) 

PORT«(POLf 1 }4P00(1>AAX2( 1A)*4X1( lA) )«S( IN) 

P|RT«fP)LU )4P11I t }*AX2( 1A>*AX1( lA) )*Yf IN) 

IXsI-l 

I.4CT(  l)))41. 

7777  PRINT  20023. I X.Z.Gl f I ) .St  IN) . Vi  I N) .POT .PI T.PORT ,PI RT 
PRINT  1757 

1757  FORMATiSSH  T  AND  XsSTATIONS  ALONG  CHORD  XHERE  X= .5( I- COSt T ) ) 

1  /SSH  GSTAR  ^CIRCULATION  DISTRIBUTION  (E02.3->S) 

1  /SSH  NCO  AND  NCIsNONLlNEAR  CORRECTION  OUTSIDE  A><0  INSIDE 

1  SSH  DUCT  RESPECTIVELY  tE03.2-3) 

1  /SSH  PO  ANO  P1«N0NL|NEAR  PRESSURE  DISTRIBUTION  OUTSIDE  AND 

1  SSH  INSIDE  THE  DUCT  RESPECTIVELY  FUR  PAr^o  DEGREES 

1  /SSH  POP  ANO  PIP  sNONLINEAR  PRESSURE  DISTRIBUTION  OUTSIDE 

1  SSH  ANO  INSIDE  THE  DUCT  RESPECTIVELY  FOR  PA= INPUT 

1) 

PRINT  103 
8889  CONTINUE 

DIMENSION  AAT(SO).XXTf 50).ZZT(50) 

READ  26000. ANT 
IF( ANT) 1801 .1802. 1801 
1801  READ  26000. ANX.ANZ 
NNX«ANX 
NNZ=ANZ 
NMT«I 

AATi 1 )3l.O 

READ  26000.(XXT( I).I«1.NNX) 

READ  2vS000.CZZr(l).I«l.NNZ) 

26021  FORMAT' 318) 

DO  26022  JA^l.NNT 
DO  26022  3X«t.NNX 
PRINT  ITSO.XXTiJX) 

DO  26024  JZ^I.nNZ 
ERASE  VACeVRG.VAO.VRQ 
DO  26023  Isl.l81.S 


ais2.*izzT( jz>-AAr( JA) }-t*>crt i ) 

B2=H«»Z«0I«»2+CXXTC JX)-1«>««2 
B5=M»*2»01*A2*’JXXT(  JX)*l«»*»2 
'’n-SORTF(XXT<  JX>  ) 

^4=B3««3 

AKsSQRTF(4.*XXT( JX)/CS) 

IF<AK-}.t  26032.26030*26030 
26032  IFIAK'.o?)  2603u. 26039.26039 

26038  CALL  ELL I P( X 1 . AK . X3. X4 .BF .  BS) 

GO  TO  26037 

26039  AK2=AK««2 

CALL  YSELL( AK2.BS.BF ) 

CO  TO  26037 
26030  8S=I.O 

26037  VAG=VA6v(-H/(6.2832«B3»»*6I ( I »«AK»( BF-BS'2. ♦( XXT( JX)- I . ) «BSA62 » ACT 
tH( I )«SMt I >AOX 

VRC=VMG4CH/(6.2832«B4} )«GU I ) aAKAHAB I • (BF-BS-2* aAKTI JX ) *85/82) ACTri 
1C  I )*SH( I )*OX 

VAQ=VA04(-H/C3.14t6*B3))*SU(l )*AK*2.*H*0t«BS*SMf 1 )*DX/B2 

26023  VRO=YR04(-H/(3.14t6*B4))*SUC  t  > AAKtC BF-BS42.*XXT ( JX ) • C XX T ( JX )-l • ) *6 
1S/82)*SM( 1) AOX 

VASVAG4VAO 

VR=VRG4VR0 

26024  PRINT  26009 >2ZT( JZ ) .VA.VR 
26022  CONTINUE 

1758  FORMAT (SSHl  INOUCEO  VELOCITY  INSIDE  THE  DUCT  AT  A  DUCT  RADIUS  OF 

ie9.4///SSM  Z  VA  VR  ) 

26000  FORMAT (9F8. 6) 

26009  FORMAT! 7E IS. 6) 

20024  F0RMAT.19I8> 

20025  FORMAT C9Fe. 6) 

801  F0RMATt//8ElS.6) 

103  FORMAT! IKl) 

29023  FORMAT! 18.8E13.4) 

;802  CALL  END  JOB 
RETURN 
END 


FOR 

SUeftOUTINE  VSELLI  EKS  «  ELE  ,  ELK) 

MS  *  1.0  •  EKS 

ALA  -  1.38629  >  fO.S  •  LOGFIPKS)) 

ELK  ■  ALA  ♦  ((ALA  -  1.01  •  0.25  •  FKS) 
ELE  •  1.0  ♦  (0.5  •  (ALA  ~  0.5MMS) 
RETURN 
END 
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